Detailed measurement of the e^e pair continuum in p + p and Au + Au collisions at 
^SArAr=200 GeV and implications for direct photon production 
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PHENIX has measured the e'*'e~ pair continuum in y^sjvjv=200 GeV An + Au and p+p colhsions 
over a wide range of mass and transverse momenta. The e^e~ yield is compared to the expectations 
from hadronic sources, based on PHENIX measurements. In the intermediate-mass region, between 
the masses of the (j) and the J/tp meson, the yield is consistent with expectations from correlated 
cc production, although other mechanisms are not ruled out. In the low mass region, below the 
(j), the p + p inclusive mass spectrum is well described by known contributions from light meson 
decays. In contrast, the Au + Au minimum bias inclusive mass spectrum in this region shows an 
enhancement by a factor of 4.7 ± 0.4"*^* ± 1.5'^"* ± 0.9™°'*'='. At low mass (mee < 0.3 GeV/c^) 
and high pr (1 < Pt <5 GeV/c) an enhanced e"'"e~ pair yield is observed that is consistent with 
production of virtual direct photons. This excess is used to infer the yield of real direct photons. In 
central Au + Au collisions, the excess of the direct photon yield over the p + p is exponential in pr, 
with inverse slope T = 221 ± 19''''" ± 19'"^"* MeV. Hydrodynamical models with initial temperatures 
ranging from Tinit — 300-600 MeV at times of 0.6-0.15 fm/c after the collision are in qualitative 
agreement with the direct photon data in Au + Au. For low pT < 1 GeV/c the low-mass region 
shows a further significant enhancement that increases with centrality and has an inverse slope of 
T ~ 100 MeV. Theoretical models under predict the low-mass, low-pr enhancement. 



PACS numbers: 25.75.Dw 

I. INTRODUCTION 

Experimental results from the Relativistic Heavy Ion 
Collider (RHIC) have established that in Au + Au col- 
lisions at y^SNN = 200 GeV matter is created with very 
high energy density 0| , as indicated by the large energy 
produced transverse to the beam direction as well 
as by the large energy loss of light ^ and heavy 
quarks ^, and is thermalized rapidly, as indicated by 
the large elhptic flow of these partons j|, 0, |[ ^- Such 
a high density thermalized medium is expected to emit 
thermal radiation in the form of direct photons and 
dileptons. 

Electron-positron pairs, or dileptons in general, are 
excellent tools for studying collisions of heavy ions at 
ultra-relativistic energies. Since they are not affected by 
the strong interaction, and therefore can escape from the 
dense medium without final state interaction, dilcpton 
spectra can probe the whole time evolution and dynam- 
ics of the collision. Dileptons can also be used to study 
the properties of low- mass vector mesons p, cu, and 4> in 
the medium, since their lifetime is shorter or similar (0) 
to that of the medium. Their mass and width inside the 
dense medium can be directly measured through their 
dilepton decay channels, and thereby one can study the 
effect of chiral symmetry restoration on these mesons. 
Furthermore, production of photons can be measured 
through their conversion to dileptons. 



As schematically shown in Fig. the dilepton spec- 
tra can be classified into the high-mass region (HMR; 
m > 3.2 GeV/c^) from J/ijj mass and above, the low- 
mass region (LMR; m < 1.2 GeV/c^) from the cj) mass 
and below, which is further subdivided in LMR I and 
LMR II as described below, and the intermediate mass 
region (IMR; 1.2 < m < 2.9 GeV/c^) between them. 

In the HMR hard scattering on partons in colliding nu- 
clei produces dileptons through the DrcU-Yan pro- 
cess {qq ~* Z+Z^) and correlated semi-leptonic de- 
cays of heavy quark pairs {bb ,cc. 
Dileptons from these hard processes are expected to 
dominate in the HMR since their mass spectra are 
harder than that from other possible sources. Thus 
dileptons in the HMR probe the initial stage of the 
collision ||ll[. Charmonia {J/ipjip') and Upsilons 
arc in this mass region and deconfinement jl^ and 
recombination Q effects can be studied from 
their yields. Little contribution from thermal radia- 
tion is expected in the HMR at RHIC energies [0 . 

In the IMR theoretical models predict that dileptons 
from the thermalized deconfincd phase, the quark 
giuon plasma (QGP), are the dominant source of 
dileptons ^ The measurement of ther- 
mal dileptons from QGP can be used to determine 
the initial temperature of the matter. Here a com- 
peting source of dileptons is the semi-leptonic decay 
of c and c, correlated through flavor conservation. 
The continuum yield in this mass region is sensitive 
to the energy loss of charm quarks in the medium. 
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In the LMR dilepton production is expected to be 
dominated by in-medium decay of p mesons in 
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FIG. 1: (Color online) Dilepton spectrum as a function of 
mass and transverse momentum from a simulation of hadron 
decays. The high-mass region (HMR; rriee > 3.2 GeV /(?) goes 
from J/tj} mass and above, the low-mass region (LMR; rriee < 
1.2 GeV/c^) from the (f> mass and below, and the intermediate 
mass region (IMR; 1.2 < m^e < 2.9 GeV/c^) between them. 
In the LMR at low-pT (II), dilepton production is expected to 
be dominated by the hadronic gas phase. Part of the LMR, 
where pT ^ mu, specifically rriee < 0.3 GeV/c^ and pT > 
1 GeV/c (I), is the quasi-real virtual photon region. The z axis 
shows the dilepton yield from the hadron decays according to 
the color scheme plotted on the right. 



the hadronic gas phase. |l^, |T^, The p has 

a strong coupling to the tttt channel, and its hfe- 
time (1.3 fm/c) is much shorter than the expected 
hfetime of the hadronic gas. The shape and the 
yield of the mass spectrum can test predicted in- 
medium modifications of the properties (the mass 
and the width) of p mesons due to chiral symmetry 
restoration ||2^ . Dileptons can also arise from other 
hadronic sources. These dilepton sources compete 
with a large contribution of e+e~ pairs from Dalitz 
decays of pseudoscalar mesons {Tr^,ri,r]') and de- 
cays of vector mesons (p, uj, (p). 

In the LMR I (marked with 1 in Fig. P is the quasi- 
real virtual photon region, where the pr of the 
dilepton is much greater than its mass {pr ^ rnu). 
Any source of real photons must also emit vir- 
tual photons which convert to low-mass e+e~ pairs. 
These low-mass pairs are produced by a higher or- 
der QED correction to the real photon emission 
process, and their yield is related to that of real 
photons. Thus e^e^ pairs in this region provide an 
alternative method for measuring direct photons. 
The measurement of the direct photon yield using 
low-mass lepton pairs was first used at the CERN 
ISR |2^. UAl observed that the low-mass dimuon 
cross section was consistent with the so-called "in- 
ternal conversion" of direct photons p3]. 



The discovery of a large enhancement of the dilep- 
ton yield in the LMR in ion-ion collisions by HE- 
LIOS/3 H and CERES (2| at the CERN SPS has trig- 
gered a broad theoretical investigation of modifications 
of properties of hadrons in a dense medium and of how 
these modifications relate to chiral symmetry restoration 
pl| , 26, |2^. These studies advanced with the availabil- 
ity of more precise data from NA60 and CERES ||, ||| 
and HADES Most theoretical studies suggest that 
in-medium modifications of the p meson, with its short 
lifetime and its strong coupling to the tttt channel, are 
primarily responsible for the enhancement. 

In the LMR II (marked with 11 in Fig. |l|) the CERES 
data show that the enhancement increases signif- 
icantly faster than linearly with charged particle 
density and is concentrated at very low pair-p^ [ p9| ■ 
This behavior is consistent with the interpretation 
that the excess is due to annihilation processes. 
NA60 recently confirmed the excess of dileptons in 
the LMR in In-|-ln collisions at 158 A • GeV with 
a high statistics dimuon measurement NA60 
also observed that the inverse slope parameter Tcff 
of the pair-px spectra rises with dimuon mass in 
the LMR [Mi. 



An enhanced yield was also observed in the IMR 
by HELIOS/3 ||, ||, NA38/50 [||, ||, and NA60 
p8| , |3l[ |. NA60 data suggest that the enhancement can- 
not be attributed to decays of D mesons but may re- 
sult from prompt production, as expected for thermal 
radiation Furthermore NA60 measures the inverse 

slope parameter of pair-p^ spectra in the IMR around 
190 MeV, independent of mass and significantly lower 
than those found at masses below 1 GeV/c^ ^ ||. 

The PHENIX experiment at the Relativistic Heavy Ion 
Collider (RHIC) has measured the dilepton continuum 
in a new energy regime, ^JsJTn — 200 GeV, for p + p 
and Au -t- Au collisions. In this article we present results 
from Au -I- Au collisions taken in 2004 and p+p collisions 
taken in 2005. We show the results in the LMR and in the 
IMR as well as the result of direct photon measurement 
from the analysis of quasi-real virtual photons. The main 
physics results in the LMR and IMR in Au + Au and p+p 
have been reported in |^ and , respectively, and the 
results of the virtual photon analysis have been reported 
in . New results on the centrality and pT dependence 
of the e^e~ pairs in the LMR are presented in this paper. 

This article is organized as follows. Sec. |l^ describes 
the PHENIX detector system related to the analysis. 
Sec. Ill presents the analysis details including the system- 
atic uncertainties. Sec. [V describes the methods used to 



calculate the pair yield expected from hadronic decays. 
Sec. shows the e+e" results as a function of mee and 
Pt, which are then discussed in Sec. ^ and c omp ared 
to available theoretical predictions. Finally Sec. VH pro- 
vides a conclusion. 
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II. PHENIX DETECTOR 

A detailed description of the complete PHENIX detec- 
tor system can be found elsewhere [|9[ |l], ||, ||, 
Here we describe the parts of the detector system 
that are used in this analysis, namely, two global detec- 
tors and two central arm spectrometers. The global de- 
tectors are the beam-beam counters (BBC) and the zero- 
degree calorimeters (ZDC). Each central arm covers the 
pseudorapidity range \r]\ < 0.35 and an azimuthal angle 
of 7r/2, and includes a drift chamber (DC) and multi-wire 
proportional pad chambers (PC) for the charged particle 
tracking, a ring-imaging Cerenkov counter (RICH) for 
electron identification, and an electromagnetic calorime- 
ter (EMCal) for energy measurement and further elec- 
tron identification. Figure ^ shows a beam view of the 
PHENIX detector. 




Beam View East 



FIG. 2: (Color online) Beam view (at 2 = 0) of the PHENIX 
central arm detector in Run-4 Au + Au and Run-5 p+p- The 
detectors used in the present analysis are the drift chamber 
(DC) and the multi-wire proportional pad chamber (PC) for 
charged particle tracking, the ring-imaging Cerenkov counter 
(RICH) for electron identification, and the electromagnetic 
calorimeter (lead-scintillator (PbSc) and lead-glass (PbCl)) 
for energy measurement. 



A. Global Detectors 

The BBC and the ZDC measure the start time and the 
collision vertex position z-vertex along the beam axis and 
are used to determine the centrality of the collision [^ . 
They also provide first level trigger information. 

The BBC consists of two sets of 64 Cerenkov counter 
modules, located ±1.44 m from the nominal interaction 
point along the beam line and measure the number of 
charged particles in the pseudorapidity region 3.1 < |ry| < 
3.9. They provide the start-time of the collision with a 
resolution of 20 ps, which gives the ^-vertex position with 



a resolution of ~ 2 cm in p + p collisions. For Au -I- Au 
central collisions we achieve a resolution of ~ 0.6 cm. 

The ZDC consists of two hadronic calorimeters, located 
'-^ 18 m from the interaction point, that measure neutral 
energy emitted within ^ 2 milliradians of the beam direc- 
tion. This energy is carried by neutrons produced either 
by Coulomb dissociation of the beam particle or by evap- 
oration from beam spectators. The energy resolution of 
the ZDC is 6E/E ~ 218%/ ( GeV) Q. 

The centrality of Au -I- Au collisions is determined by 
the correlation between the BBC charge sum and the 
ZDC total energy @. 

B. Central Magnet 

The Central Magnet (CM) is an axial field magnet 
energized by two pairs of concentric coils, roughly in a 
Helmholtz configuration, which can be run separately, 
together, or in opposition so that the momentum and 
the charge of a particle can be determined by its bend- 
ing curvature | ^l[ |. In the mode in which both coils are 
running in the same direction, the single particle mo- 
mentum resolution is better than 1% between 0.2 and 
1 GeV/c. During the Au -I- Au measurement in Run-4 
and the p+p measurement in Run-5, the field component 
parallel to the beam axis had an approximately Gaussian 
dependence on the radial distance from the beam axis, 
dropping from 0.9 T at the center to 0.096 T (0.048 T) 
at the inner (outer) radius of the DCs. The total field 
integral is J B ■ dl ^ 1.15 Tm. 

C. Tracking Detectors 

The drift chambers (DCs) and the pad chambers 
(PCs) in the central arms measure charged parti- 
cle trajectories in the azimuthal direction to determine 
the transverse momentum (pr) of each particle. The DC 
provides the most precise measurement of particle tra- 
jectories in the plane perpendicular to the collision axis. 
The first layer of PC provides the most precise measure- 
ment of the track space point along the collision axis. 
Additional layers principally supply pattern recognition 
support. 

The DCs are located in the radial region 2.02 < r < 
2.46 m. Each DC volume consists of 20 sectors, each of 
which covers 4.5 degrees in azimuth and |r?| < 0.35. Each 
sector has six types of wire modules stacked radially: XI, 
Ul, VI, X2, U2 and V2. Each module is further divided 
into 4 drift cells in the (j) direction. A plane of sense 
wires is at the center of each drift cell, with 2 to 2.5 
cm drift space on either side. The X wires run parallel 
to the beam axis and measure the particle trajectory in 
the r-(j) plane. The U and V wires have a stereo angle 
of about 6.0 degrees relative to the X wires in order to 
measure the z-coordinate of the track. Each X- and U, 
V-stereo cell contains 12 and 4 sense wires, respectively. 
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The single X wire position resolution is 150 /xm. The 
intrinsic tracking efficiency of the X modules is greater 
than 99 %. 

The pad chambers (PC) are multi-wire proportional 
chambers that form three separate layers. They deter- 
mine space points along the straight line particle tra- 
jectories outside the magnetic field. The first PC layer 
(PCI) is located between the DC and the RICH, the sec- 
ond layer (PC2) is placed behind RICH (west arm only) 
and the third layer (PCS) is located in front of the EM- 
Cal. PCI and the DC, along with the z— vertex position 
measured by the BBC, are used in the global track re- 
construction to determine the polar angle of each charged 
track. The position resolution is ±1.7 mm for PCI along 
the wire (z-direction). 

Helium bags were installed between the beam pipe and 
the DCs to reduce the conversion material prior to the 
first tracking layer to ~ 0.4% of a radiation length. The 
material budget is known with an uncertainty of ~ 5%. 



The PbSc is a Shashlik-type sampling calorimeter 
made of alternating tiles of lead and scintillator. It con- 
sists of 10.5 X 10.5 X 37 cm^ (18.2 Xq) modules, con- 
structed of alternating layers of 1.5 mm thick lead, re- 
fiecting paper, and 4 mm thick polystyrene-based scin- 
tillator. Each module is divided into four equal towers, 
from which the light is collected separately by scintillat- 
ing fibers. Each PbSc sector consists of 36((/)) x 72 (z) = 
2592 towers. The nominal energy resolution is 6E/E ^ 
4.5% e 8.3%/y^E ( GeV). 

The PbGl is a Cerenkov counter that measures the 
light emitted by the particles in an electromagnetic 
shower and collected by one PMT at the back end. 
Each PbGl sector consists of 4608 4.0 x 4.0 x 40.0 cm^ 
(16 Xq) modules made of lead-glass crystals. The PbGl 
has a nominal energy resolution of SE/E 4.3% ® 



7.7%/ y/E ( GeV). 



III. ANALYSIS 



D. Ring-imaging Cerenkov Counter (RICH) 

The RICH is a threshold gas Cerenkov detector 
and is the primary detector to identify electrons in 
PHENIX It is located in the radial region 2.5 < r < 
4.1 m, just outside PCI. Each arm contains spherical 
mirror panels (0.53% of a radiation length) which focus 
Cerenkov light onto two arrays of 80{(j)) x 16(z) = 1280 
PMTs. The PMTs are located outside the acceptance, 
on either side of the RICH entrance window, and are 
shielded to allow operation in a magnetic field up to 100 
Gauss. The Cerenkov radiator gas, CO2 at atmospheric 
pressure, has n = 1.000410 (7 ~ 35) which corresponds 
to a momentum threshold of 20 MeV/c for an electron 
and 4.65 GeV/c for a pion. 

The average number of hit PMTs per electron track 
is about 5, and the average number of photo-electrons 
detected is about 10. Simulation studies show that pion 
rejection by the RICH alone, for isolated tracks, is limited 
by the production rate of "co-linear" delta electrons to 
one part in 10^. However in high multiplicity collisions 
pion tracks may be mistaken for electrons via overlap of 
their trajectory with a true electron's ring. This effect 
worsens the RICH-alone pion rejection to roughly one 
part in 10'^ for central Au -I- Au collisions and requires 
additional cuts in the offline analysis as described below. 



E. Electromagnetic Calorimeter (EMCal) 

The EMCal |Q provides a measurement of the ener- 
gies and the spatial positions of photons and electrons. 
Each arm consists of four rectangular sectors in ip: the 
two bottom sectors of the east arm are lead-glass (PbGl) 
calorimeters; whereas, the remaining sectors are lead- 
scintillator (PbSc) calorimeters. The radial distance from 
the z axis is 510 cm for the PbSc and 550 cm for the PbGl. 



In this Section we present all steps of the data anal- 
ysis. We start by introducing the data set, the event 
selection for p + p and Au + Au and the centrality def- 
inition for Au + Au collisions (III A). We present the 
single electron analysis, including_ track reconstruction 



(IIIB) and electron identificat ion ( IIIC). We present the 
details of the pair an alysis ( illD ) including pair cuts 
and p hoton r ejection (HID 1), combinatorial and corre- 
lated (HIDJ) background subtraction. An alternative 
method to subtract combinatorial and correlated back- 
Section 



ground together is described in HID 3 



HIE 



[II F 



presents the final raw mass spectrum. In Section 
we present the mass spectrum obtained from the anal- 
ysis of runs with increased conversion material, a tech- 
nique employed to estimate the systematic uncertainty 
on the ba ckgro und subtraction. Next we describe the 



efficiency (HIG) and acceptance (IHH) corrections, trig- 
ger efficiency (Flll|) (for p + p collision s) an d occupancy 
correction (for Au -I- Au collisions) (IIIJ). Finally we 
describe the calculation of the associated systematic un- 
certainties involved in the analysis (IIIK). 



A. Data Sets and Event Selection 

The data ior p + p collisions at y/s = 200 GeV were 
collected during the polarized p + p rim in 2005. For this 
analysis two data sets were used: a reference sample of 
events selected with a minimum bias interaction trigger 
(Min. Bias) and a data set selected with a single electron 
trigger (ERT: EMCal and RICH Trigger). The Min. Bias 
trigger for p + p requires at least one hit in both the 
North and South BBC detectors in coincidence with the 
beam bunch crossing and a z-vertex position (determined 
online by the BBCs) within 38 cm: 

Min. Bias=(BBCN> 1) x (BBCS>1) x (|z| < 38 cm).(l) 
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The Min. Bias trigger cross section is measured to be 
ctbbc = 23.0±2.2 mb or 54.5 ± 5% of the total in- 
elastic p + p cross section at this center of mass energy 
crP+P = 42 ± 3 mb ||||. Data collected without requir- 
ing the BBC trigger show that the BBC fires on 79% of 
events with particles in the central arm acceptance. We 
assume that regardless of the electron px or electron pair 
Pt and invariant mass, the BBC always fired with the 
same probability of 79%. Therefore, in the p + p data 
the yield is divided by 0.79/0.545 to account for the frac- 
tion of tracks (0.79) and inelastic p + p collisions (0.545) 
missed by the Min. Bias trigger. 

The ERT trigger requires a minimum energy deposit 
of 0.4 GeV in a tile of 2x2 EMCal towers matched to 
a hit in the RICH, in coincidence with the Min. Bias 
trigger. In the active area the ERT trigger has a very high 
efficiency for electrons; it reaches approximately 50% at 
~ 0.5 GeV/c and saturates at ^ 1 GeV/c close to 100% 
(for the EMCal) and close to 90% (for the RICH). 

After applying a z- vertex cut \z\ < 25 cm, and discard- 
ing any run with unusual beam or detector conditions, 
the total integrated luminosities were 65.6 nb~^ and 2.49 
pb~^ for the Min. Bias and ERT data sets, respectively. 

The data for Au -I- Au collisions at y^ijvw = 200 GeV 
were collected during the run in 2004. Collisions were 
triggered using beam-beam counters (BBC). The Min. 
Bias trigger requires at least 2 hits in each of the BBCs 
and \z\ < 38 cm: 

Min.Bias ee (BBCN > 2) x (BBCS > 2) x (|z| < 38 ci^) 

The offiine Min. Bias trigger also requires one hit in one 
of the ZDCs. The same z-vertex cut \z\ < 25 cm as in the 
p+p data is applied offline. This corresponds to 92I3 q% 
of the Au -I- Au inelastic cross section. 

The centrality is determined for each Au -I- Au colli- 
sion by the correlation in the measurement of the BBC 
charge and ZDC energy |4^. Using simulations based 
on a Glauber model calculation j^] the average number 
of participants iVpait and the average number of binary 
collisions iVcoii associated with each centrality bin are 
determined. Table | summarizes the average iVpart and 
TVcoU and the corresponding systematic uncertainties in 
each centrality class used in the analysis. 

We analyzed a sample of 8x 10® minimum bias events, 
divided into five centrality classes (0-10%, 10-20%, 20- 
40%, 40-60%, and 60-92%) for which the number of 
events is summarized in Table m. 



B. Track Reconstruction 

The PHENIX tracking system reconstructs charged 
particles with momentum above 0.2 GeV/c with a mo- 
mentum resolution of ap^/pr — 0.7% ® 1% ■ px for pT 
in GeV/c. A track is reconstructed by 2 sets of at least 
4 hits in the XI and X2 plane separated by 20 cm in 
radial direction, i.e. in the main bend plane of the cen- 
tral magnet, using a Hough transform performed over all 



possible hit combinations. The UVl and UV2 wires pro- 
vide up to 6 measurements in the z direction, which are 
associated with the three-dimensional space point pro- 
vided by PCI. After the pattern recognition and track 
reconstruction by the Hough transform technique, the 
initial momentum vector of the track at the z-vertex is 
calculated. Each reconstructed track is then associated 
with hit information from the outer detectors (PC2, PC3, 
RICH, and EMCal). 

The transverse momentum {pt) is determined by mea- 
suring the angle a between the reconstructed particle tra- 
jectory and a line that connects the z-vertex point to 
the particle trajectory at a reference radius i?=220 cm. 
The angle a is approximately proportional to charge/pT- 
Note that this procedure assumes tracks originate from 
the vertex. As a result, tracks which originate off vertex 
are reconstructed with an incorrect momentum. Conver- 
sion pairs are reconstructed with invariant mass TOge > 
and contaminate the spectrum up to rngg ^0.3 GeV/c^ 
(see Section pIDH ). 

Because charged particles are defiected in the az- 
imuthal direction by the magnetic field, the single-track 
acceptance depends on the momentum and charge of the 
particle, and also on the radial location of the detector 
component (DC and RICH). The acceptance for a track 
with charge q, transverse momentum px and azimuthal 
emission angle (j) can be described by the logical AND of 
these conditions: 



(j^-miYi <4> + q-^ < (/'max 
PT 



, , , fcRic: 

Pmin <<P + q 

PT 



(3) 



where fc^c and fci^icH represent the effective azimuthal 
bend to DC and RICH (fcoc = 0.206 rad • GeV/c and 
fcRiCH = 0.309 rad • GeV/c). One arm covers the region 



from 



11- 
16' 



'-n to 



rvr, the other arm from 



16 • 



Results in Sec. [V] will show the dilepton invariant mass 
spectrum "in the PHENIX acceptance," where the data 
will be compared to the expectations filtered according 
to this simple parameterization of the acceptance. 



C. Electron Identification 

Electrons in the range 0.2 < < 20 GeV/c are iden- 
tified by hits in the Ring Imaging Cerenkov detector 
(RICH) and by matching the momentum with the en- 
ergy measured in the electromagnetic calorimeter (EM- 
Cal) Q . Specifically we consider the following variables 
for electron identification (elD), as summarized in Table 

0- 

track quality A bit pattern representing the reconstruc- 
tion quality of the track. If the track is recon- 
structed by both of the XI and X2 sections of the 
DC and is uniquely associated with hits in U or 
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TABLE I: Nr^it, NcoU for Au + Au collisions at -^sjvjv ~ 200 GeV with the corresponding uncertainties derived from a Glauber 
calculation pi and the number of events and pairs for each centrality class. Note that the uncertainties are correlated P8|. 



Centrality class 




\ ^Vcoii/ (systj 


J 'events 


Signal Pairs 


U-iU 70 


■ilo.l [6.6) 


900.4 (93. 6 j 


O.D X iU 


9.Z X iO 


10-20 % 


234.6 (4.7) 


602.6 (59.3) 


8.6 X 10'^ 


6.6 X 10* 


20-40 % 


140.4 (4.9) 


296.8 (31.1) 


1.7 X 10* 


8.1 X 10* 


40-60 % 


59.95 (3.6) 


90.70 (11.8) 


1.7 X 10* 


3.3 X 10* 


60-92 % 


14.50 (2.5) 


14.50 (4.00) 


2.9 X 10* 


1.1 X 10* 


0-92 % 


109.1 (4.1) 


257.8 (25.4) 


8.1 X 10* 


28.3 X 10* 


p + p (Min. Bias) 


2 


1 


1.5 X 10^ 


1.4 X 10* 


p + p (ERT) 


2 


1 


2.7 X 10* 


22.8 X 10* 



V stereo wires, the value of quality is 63 (in case 
a unique PCI hit is found too) or 31 (in case the 
PCI hit is found but ambiguous). If there are no 
UV hits found, but a PCI hit is, quality is 51. 

EMCal match (ca^) Displacement in between the po- 
sition of the associated EMCal cluster and the pro- 
jection of the track onto the EMCal. The quantity 
is measured in units of momentum-dependent res- 
olution. For example, A0 < 2 means that the posi- 
tion of the associated EMCal cluster in is within 
2(7 of the projected track position. The particle hit 
position of an EMCal cluster is particle-species de- 
pendent due to different shower shapes. Here the 
parameterization has been optimized for electrons. 

EMCal match (ctaz) Analogous to the previous vari- 
able, for the z coordinate. 

nO Number of hit RICH PMTs in an annular region with 
an inner radius of 3.4 cm and outer radius of 8.4 cm 
around the track projection on the RICH. The ex- 
pected radius of a Cerenkov ring emitted by an elec- 
tron is 5.9 cm. 

chi2/npe0 A x^-like shape variable of the RICH ring as- 
sociated with the track divided by the number of 
photo-electrons measured in a given ring (npeO). 

RICH match The displacement of the RICH ring center 
from the projected track position. Units are cm. 

E/p or dep A variable quantifying energy-momentum 
matching. This variable is calculated as dep — 
{E/p — l)/aE/p, where E is the energy measured 
by EMCal; p is the momentum of the track; and 
Ue /p is the standard deviation of the Gaussian- like 
E/p distribution calculated for electrons. 



TABLE II: Electron ID cuts used in the Au -I- Au and p + p 

analyses. 



elD cuts 



for Au -I- Au 



for p + p 



track quality = 



63 



A4> ' "Az 

nO > 
chi2/npe0 < 
RICH match < 
dep > 

# > 



I 31 II 
3.0 

2 
10.0 
5.0 
-2.0 



51 



63 



31 II 
5.0 

1 

15.0 
10 

0.5 



51 



10=" 



10" E 



10'' E 



min. bias Au+Au \/s^ = 200 GeV 

• All charged tracks 

o + RICH 
— + random RICH 




I I I I I I I I I I 

0.2 0.4 0.6 



FIG. 3: (Color online) E/p distribution in minimum bias 
Au + Au for all charged tracks and for tracks after applying 
all the RICH cuts in Table ^ except the E/p or dep. The 
contribution from randomly associated hadrons is shown by 
the filled histogram. 
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Figure ^ shows the E/p distribution for all charged 
tracks and for electron candidates, i.e. tracks which ful- 
fill all the RICH elD cuts except the E/p or dep. While 
the distribution of all charged tracks shows no clear 
electron peak, requiring the elD cuts greatly improves 
the signal-to-background ratio. However, there still re- 
mains some background underneath the peak even below 
Pt < 4.9 GeV/c, the Cerenkov threshold for pions. This 
background, due to random coincidences between hadron 
tracks and hits in the RICH, is estimated by swapping the 
north and south sides of the RICH in software, and re- 
constructing the track matching to the RICH once again. 
This contamination ranges from 2% in p + p to 30% in 
the most central Au -|- Au collisions. It is ~24% for 
Min. Bias Au -|- Au collisions. This pion contamina- 
tion in Au + Au con tribute s to the large combinatorial 
background (see Sec. HI D S ). 



D. Backgrounds 

The source of any particular electron or positron 
in an event is unknown; therefore, all electrons and 
positrons are combined into a foreground of pairs, like- 
sign 7V_|_+,7V and unlike-sign This results in a 

large combinatorial background which must be removed. 
In the following we will use the notation for the fore- 
ground N±± = iV++ + Af__ and for the background 

B±± = -B++ -I- B The analysis steps to achieve this 

are outlined here and presented in detail in the subsec- 
tions below. We can distinguish our background sources 
in two types: 



Type I background ( [II D 1 ) consists of two classes of 
fake combinations that can be identified on a pair-by-pair 
basis: 

• Overlapping pairs are fake electron pairs that arise 
from overlapping hits in the detectors, mostly in 
the RICH. 

• Photon conversions are fully reconstructed pairs 
originating from photon conversions in the detector 
material are removed by a cut on the orientation of 
the pairs in the magnetic field. 



Type II background (HID 2) consists of all those pairs 
that cannot be identified on a pair-by-pair basis and are 
therefore removed statistically: 

• Combinatorial background B'^°™^ arises from all 
the combinations where the origin of the two elec- 
trons is totally uncorrelated. 

• Correlated background B"^"" occurs if there are two 
e~^e~ pairs in the final state of a meson or when two 



^ Data tables for this and other data plots are available at 
https: / /www.phcnix.bnl.gov/WWW/p/info/ppg/088/datatables. 



hadrons, either within the same jet or in back-to- 
back jets, decay into electron pairs. 

Since accurate background subtraction is essential for 
this analysis, we have developed two independent meth- 
ods to subtract the type II background. In the first 
method we calculate the shapes of combinatorial and cor- 
related background with event mixing or simulations and 
use the yield of the like-sign s pectra for the normaliza- 
tion. In the second method ( III D 3 ) we use measured 
the real like-sign distributions, corrected for the accep- 
tance difference, without making any assumption about 
the mass dependence nor about the decomposition of the 
background into correlated and uncorrelated component. 



1. Type I Background 

Overlapping Pairs Fake electron pairs can be cre- 
ated if two particles are in close proximity in any of the 
detectors. These correlations within an event are partic- 
ularly noticeable when two tracks share the same RICH 
ring. This issue can be illustrated with a simplified model 
of the RICH based on spherical mirror optics. In this 
case, tracks that are parallel to each other while pass- 
ing through the RICH radiator, i.e. after they have been 
bent in the magnetic field, share the same search region 
for Cerenkov light in the RICH PMT plane. Therefore 
an overlapping pair is created whenever a track after the 
field is parallel to a true electron. These overlapping 
pairs typically have a small opening angle and are there- 
fore reconstructed with small invariant mass. Because 
like-sign pairs are bent in the same direction, in contrast 
to unlike-sign pairs which are bent in opposite directions, 
like-sign overlapping pairs have smaller mass than unlike- 
sign overlapping pairs. Figure ^ shows the invariant mass 
distribution for all like- and unlike-sign pairs, the distri- 
butions for overlapping pairs, and the distributions after 
removing overlapping pairs. 

Overlapping pairs are eliminated by applying a cut on 
the physical proximity of every pair projected to every de- 
tector. The cut value on every detector is determined by 
the corresponding double hit resolution. In the RICH the 
cut is applied at 36 cm, which is roughly twice the pre- 
dicted maximum diameter (~ 16.8 cm) of a RICH ring. 
In the Pad Chambers the cut is applied at Az < 0.5 
cm and A(p < 20 mrad. In the EMCal the cut is ap- 
plied to a 3x3 tower region around the hit. These cuts 
remove a fraction of real pairs which varies from 4% in 
the most central to 2% in the most peripheral collisions, 
estimated using mixed events. These cuts remove 10% 
more real like-sign than unlike-sign pairs. The ratio of 
like- to unlike-sign pairs lost was determined with an un- 
certainty of 50% by varying the cut values chosen. While 
the efficiency of the pair cut depends on the centrality, 
the ratio of like- to unlike-sign pairs lost was found to be 
independent of centrality. 

Because these overlapping pairs are rare, whenever we 
encounter one, we remove the entire event. This results 
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in a loss of ^ 0.08% of all events. In mixed events the 
same cuts are applied and whenever an event is discarded 
due a pair cut, another event is generated with the same 
electron multiplicity and in the same centrality, z-vertex 
and reaction plane class. 




FIG. 4: (Color online) Invariant mass distribution for all like- 
and unlike-sign pairs. Overlapping pairs are shown separately. 
Also shown are the distributions after the overlapping pairs 
are removed. 

Photon Conversions Since the tracking algorithm 
assumes that all tracks originate from the collision ver- 
tex, pairs from photons that convert off- vertex are recon- 
structed with an artificial opening angle, which leads to 
an invariant mass that increases with the radius at which 
the conversion occurs. 

Conversion pairs have no intrinsic opening angle (i.e. 
their opening angle is exactly zero at the point of cre- 
ation), they are bent only in the azimuthal direction by 
the magnetic field, which is parallel to the beam axis z. 
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FIG. 5: (Color online) Comparison of the (pv angle distribu- 
tions in Monte Carlo and data in the mass range 10 < niee < 
30 MeV/c^. The vertical red line indicates the cut value used 
to eliminate photon conversions in this mass bin. 



We can define unit vectors u in the direction of the pair 
momentum and v perpendicular to the plane defined by 
the pair 



P++P- 
\P+ + P- 1 
P+ X p- 



(4) 
(5) 



where p± — p±/\p±\. We can define the orientation of 
the actual opening angle as 



W ~ U X V 



(6) 



We can also define the expected orientation of the open- 
ing angle for conversion pairs 



Wc 



(7) 



Finally we can define (f>v as the angle between these two 
vectors 



cos (pv 



(8) 



For pairs originating from photon conversions 4>v is zero. 
(By consistently ordering positive and negative tracks 
within the pair we avoid (j)v = tt as a solution for pho- 
ton conversions.) In contrast, e'^e~ pairs from hadron 
decays, as well as combinatorial pairs, have no preferred 
orientation. Figure || presents a comparison of the 4>v an- 
gle distributions of real data and Monte Carlo simulated 
data in a mass bin, 10 < mee < 30 MeV/c^, dominated 
by photon conversion in the beam-pipe. In the simula- 
tions we can distinguish which e+e~ pairs originate from 
a TT*^ Dalitz decay (dotted-dashed line) and which origi- 
nates from photon conversion (dotted line). The simula- 
tions show that the distribution for all unlike-sign pairs 
originating from photon conversion is strongly peaked at 
(j)v — 0. In contrast, simulated pairs originating from 
tt'^ Dalitz decays have no strongly preferred orientation. 
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FIG. 6: (Color online) Invariant mass spectrum of all the 
unlike-sign pairs after subtraction of combinatorial pairs in 
Min. Bias Au + Au collisions at y/s^ = 200 GeV. The filled 
histogram shows the pairs removed by the (pv angle cut. 



2. Type II Background 

After removing type I background, the unlike-sign fore- 
ground spectrum measures the physical signal plus 

background, while the like-sign spectra N+^,N mea- 
sure only background. We have developed two methods 
to measure and subtract the unlike-sign background 

• One solution is to use a mixed-event technique, 
which combines tracks from different events. With 
this method the background has much larger statis- 
tics than the foreground. The accuracy in the 
determination of the shape of the background is 
tested by comparing the like-sign distribution in 
real and mixed events. We find good agreement 
between real and mixed-events like-sign spectra in 
some regions of the (mecPr) plane, while in oth- 
ers they clearly deviate. This indicates that not 
all the type II background is of uncorrelated ori- 
gin, but ther e are al so some correlated pairs in the 
background (HID 2) that one needs to separately 
account for. 



The sum of the Monte Carlo data can be compared to 
the signal, which in this mass region contains conversion 
photons, as well as tt^'s. Figure || shows that the agree- 
ment between the signal and the simulation is good. The 
width of the (j}y peak for photon conversions increases 
with the path length of the e~^e~ pair in the magnetic 
field, where residual field in the polar direction as well 
as multiple scattering cause the pair to lose its perfect 
alignment perpendicular to the z axis. 

The contribution from conversion pairs as a function of 
the (mis-reconstructed) invariant mass is shown in Fig- 
ure H (filled histogram). As the (mis-reconstructed) mass 
is essentially proportional to the radius where the conver- 
sion happens, the mass spectrum of those pairs allows a 
"tomography" of the material in the spectrometer. The 
peaks correspond to the conversions in the beam-pipe 
material (r = 4 cm, or ruee — 20 MeV/c^) and detector 
support structures (r = 25 cm, or rriee = 125 MeV/c^). 
Conversions in the He bag generates pairs with niee ^ 
0.3 GeV/c^. For this value of mass the corresponding 
radius would be the entrance window of the DC. At this 
point though the electrons do not bend anymore because 
the region is field-free, and are therefore removed with a 
Pt cut {pt <20 GeV/c) on the single electrons. The (py 
resolution improves for increasing conversion radius be- 
cause electrons are less modified from their original direc- 
tion by multiple scattering or the residual polar field. The 
cut on (j)v (indicated in Fig. |^ by the line at (j)v > 0.25 
for rriee < 30 MeV/c^) is reduced for larger masses due to 
the improved resolution of (fiy at larger radii. By varying 
the cut values on real pairs we estimate that our 4>v cut 
removes more than 98% of the conversion pairs with a 
mass-dependent efficiency of more than 90%. The uncer- 
tainty of on the final e'^e~ pair signal is 6%. 



• Another solution is to still use the measured like- 
sign spect ra and correct them for the different ac- 
ceptance ( HI D 3 ). This solution does not require 
any assumption on the decomposition of combi- 
natorial and correlated since the like-sign spectra 
measures all the backgrounds simultaneously. In 
experiments with equal acceptance for electrons 
and positrons, the background can be measured di- 
rectly through the geometric mean of the like-sign 

pair distributions 2^yN++N With this method 

the background has similar statistics as the fore- 
ground. In PHENIX however the mass and pr de- 
pendence of the background is different for the two 
charge combinations (see Equation || and Fig. 

We have used b oth these methods, which are compared 
to assign a systematic uncertainty to 



in Section [II D 4 



the background subtraction. The first method is used for 
the final analysis. 

The structure of the type II background was studied 
using minimum bias events generated with pythia 
with the branching ratio of the 7r° Dalitz decay set to 
100% to enhance the sample of e+e" pairs per event. 
All the electrons are filtered through the PHENIX ac- 
ceptance. From PYTHIA events we make real and mixed- 
event like- and unlike-sign distributions as in the real 
data. Now we analyze like-sign pairs, that contain only 
background. The Ar/i distribution shown in Fig. |^ com- 
pares the difference of azimuthal emission angle of the 
two electrons (or positrons) in real and mixed pythia 
events. If the two shapes were identical, we would con- 
clude that all the background come from uncorrelated 
sources. However the shapes clearly deviate at At/) ~ 
and A(j) tt. This indicates the presence of combina- 
tions that arise from the same jet (Acj) ~ 0) or back-to- 
back {A(j) ~ tt) jets (see Section HI D 2| ) . In addition, 
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correlations can occur if there are two e+e" pairs in the 
final state of a meson, e.g. double Dalitz decays, Dalitz 
decays followed by a conversion of the decay photon or 
two-photon decays followed by conversion of both pho- 
tons. 
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FIG. 7: (Color online) A(/!> distribution of like-sign pairs for 
real and mixed events calculated by pythia p9[ . 

Therefore, we decompose the type II background into 
two components: a combinatorial background made of 
uncorrelated pairs and a background of correlated pairs. 
The distributions for combinatorial and correlated back- 
ground are de termined with methods that are explained 
below (HID 2) for like and unlike-sign spectra. The like- 
sign background distributions are fit to the measured like- 
sign spectra (that contain only background). The same 
fit parameters are then applied to the unlike-sign back- 
ground. The normalized unlike-sign background is finally 

subtracted from the foreground of all pairs to obtain 

the signal. 

Combinatorial Background The combinatorial 
background B'^"™^ is determined with a mixed-event tech- 
nique, which combines tracks from different events with 
similar centrality, z-vertex and reaction plane. In the 
p + p data, where we use a triggered data set, the mixed- 
event pairs are constructed from the Min. Bias data set 
requiring that at least one of the two partners has ful- 
filled the ERT trigger condition. Since the tracks are 
from different events, this technique reproduces the un- 
correlated background by definition. This technique also 
allows computation of background spectra with negligible 
statistical errors. 

We compare like-sign spectra in real and mixed-event 
data to locate a region in the (rnee,PT) plane where their 
shapes agree. In this region we normalize the combinato- 
rial background spectra B"^^^ to the measured like-sign 



pairs N±±. We define: 



A, 



/jV.fl. ^++(™ee,PT) 
In.R. B++imee,PT) 



dwLeedpT 



where N.R. is the chosen normalization region. Then we 
calculate the integral of the normalized like-sign back- 
ground over the full phase-space: 

oo 

B++ = J A+ ■ B++{mee,PT)dmeedpT 


oo 

^-=J^--B-irn^e,PT)dm..dpr (10) 



The unlike-sign background is then normalized such that 
its yield equals the geometric mean of the like-sign pairs 
2v/fi++6__: 



B+^{mee,PT) 



B+_{mee,PT)dmeedpT 



B+^{mee,PT<)ll) 



Appendix ^ shows that as long as electrons and positrons 
are produced in pairs the absolute normalization of the 
unlike-sign background is given by the geometric mean 
of the observed positive and negative like-sign pairs 
2yJ N^^N , without any further assumption about ef- 
ficiencies, acceptances or probability distribution func- 
tions for the pair. Using B++ and B instead of directly 

taking iV++ and N simply avoids counting correlated 

pairs, w hich will be measured and normalized separately 
(|IIID2| ). 

The systematic uncertainty of the normalization is 
therefore determined by the statistical accuracy of the 
measured like-sign yield in the region chosen for the nor- 
malization. 

Combinatorial Background in p + p data 

In our PYTHIA 1^ studies we found (Fig. 0) that at 
A0 ~ and A(/) '-^ tt the real events background de- 
viates from the shape of uncorrelated sources, while at 
A</) ~ 7r/2 it looks consistent with the shape from mixed 
events. Considering that m^g = 2pip2(l — cos Aw), with 
Aw being the opening angle, this condition corresponds 
to a region in the (rriecPT) plane where Wee ~ Pt- 
We define this region empirically by a set of equations: 

TOee > 0.3 GeV/c^ 
niT < 1.2 GeV/c^ 
Pt/c - 1.5mee < 0.2 GeV/c^ 
Pt/c - O.TbrUee > GeV/c^ 



(12) 



. In B. N^^{mee,PT) ^ 
A- = " „ 7 -dmeedpT 



(9) 



where tot = •\/pt+"^ transverse mass of the 

pair. This region is shown in Fig. ^ by the dashed area. 

Figure ^ shows the difference between the like-sign dis- 
tributions in real and mixed-events, as a function of rrige 
and Pt- The background _B™^* is normalized to the 
foreground N±± in the normalization region (from Equa- 
tion |l^) . The absolute normalization of the unlike-sign 
combinatorial background is determined with an uncer- 
tainty given by the statistical error of the measured like- 
sign spectra in this region of 3%. The difference between 
real and mixed N±± — B'^™^ is divided by its standard 
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deviation. Figure ^ shows that in this region the back- 
ground does not deviate from the foreground by more 
than 2 X (7. The stabihty of the results has been checked 
by varying the normahzation region, and the difference 
is included in the systematic uncertainty. 
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FIG. 8: Difference between real and mixed-events like- 
sign distributions divided by its standard deviation (A'^±± — 
S±±'')/<^(iv±±-B5f>4"b). The background B'^^^ is normalized 
to the foreground N±± in the normahzation region shown by 
the dashed area. 



Combinatorial Background in Au + Au data 

Figure || shows a comparison between the like-sign dis- 
tribution from real and mixed events. The compar- 
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FIG. 9: (Color online) (a) like-sign distribution for real 
iV±± and mixed events B±±''. (b) and (c) Ratio of (iV±± - 
^conib)/^comb_ with two different scales. 
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FIG. 10: (left) Maximum cumulative fractional distance 
d of the two like-sign distributions. The distance is cal- 
culated from a lower endpoint miow to infinity. The test 
gives a maximum deviation of 0.1%. (right) Corresponding 
Kolmogorov-Smirnov p-value as a function of the lower end- 
point miow , i.e. for rries > miow The p- value increases for 
miow > 0.55 GeV/c^, reaching values of ~ 90%, confirming 
the hypothesis of compatibility of the two distributions in the 
region chosen for the normalization. 



ison shows that the mixing technique reproduces the 
mass dependence within the statistical accuracy of the 
data not only in the normalization region for the p+p 
data, but also for all masses above 0.55 GeV/c^. In 
general the larger combinatorial background produced 
in the Au -|- Au environment reduces the capability to 
distinguish between different background shapes. Also 
the agreement at high mass, which was not observed in 
p + p (see Fig. can be qualitatively explained by the 
suppression of away-side jets observed in Au + Au | |5C| ]. 

A small signal from c orrelate d background remains at 
low-masses (see Section p^II D 2 ). To quantitatively com- 
pare the mass dependence of the data to the mixed events 
we calculate the ratio (iV±± - B^l'^)/B^f' shown in 
the bottom panels of Fig. ^ and fit it with a con- 
stant above the 77 mass (0.55 GeV/c^). The result is 
(-2.59 ± 6.33) X 10-"* with xVNDF = 27.6/19. Panel 
(b) shows the entire mass range and allows to distinguish 
the signal from correlated background at very low mass. 
Panel (c) shows a zoom in the region where we fit. 

Two statistical tests are performed to test the hypoth- 



esis that the two distributions {N±± and ) repre- 

sent identical distributions: the Pearson test and the 
Kolmogorov-Smirnov test. 

The Pearson test statistic resulting from the com- 
parison of real and mixed events for ruee > 0.55 GeV/c^ 
returns a p- value greater than 0.83. 

The Kolmogorov-Smirnov test is used to decide if a 
sample comes from a population with a specific distribu- 
tion, either comparing one fluctuating distribution (the 
test) to a truth hypothesis (reference), or querying two 
fluctuating distributions as to whether they have a com- 
mon truth origin. This latter one was done in our case 
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TABLE III: Fit parameters for the mass dependence comparison of real and mixed-events like-sign pairs for different centrality 
and Pt bins. The second column reports the results of (A^±± — B§?™'')/_B§°™'' distributions to a constant and the third is 
value divided by the number of degrees of freedom. The fourth and fifth columns report the result of a statistical test 
and the corresponding p-value for rUee > 0.55 GeV/c^. The last one gives the maximum deviation of the N±± and B'^™^ 
distribution in a Kolmogorov-Smirnov test. 



Centrality 


Po 


xVndf 


X^ test 


p-value 


max dev. 


0-10% 


6.3 ± 8.8 X 10 


30.2/19 


1.05 


0.25 


0.0014 


10-20% 


-9.4 ± 1.4 X 10"" 


18.6/19 


0.97 


0.61 


0.0018 


20-40% 


-2.4 ± 1.8 X 10"^ 


18.7/19 


1.02 


0.40 


0.0034 


40-60% 


-8.5 ±4.9 X 10~^ 


21.9/19 


1.65 


0.02 


0.0071 


60-92% 


-1.8 ± 1.6 X 10~2 


21.5/14 


1.51 


0.04 


0.0321 


00-92% 


2.6 ±6.3 X 10"" 


27.6/19 


0.92 


0.83 


0.0010 


Pt < 1 GeV/c 


9.2 ±5.1 X 10"" 


18.9/18 


0.95 


0.73 


0.0011 


1< PT < 2 GeV/c 


-3.4 ± 1.6 X 10"^ 


27.9/18 


0.91 


0.84 


0.0029 


Pt > 2 GeV/c 


-9.6 ±5.4 X 10"^ 


15.2/18 


0.97 


0.63 


0.0038 



since we do not have a truth reference distribution with- 
out fluctuations. The test is based on the maximum cu- 
mulative difference between the distribution under test 
and a specified reference distribution. The left panel 
of Fig. |o| shows the maximum cumulative difference d 
of the two like-sign distributions 7V±± and 3"^^"° for 
rriee > JTiiow, where miow is the lower endpoint chosen 
between 0.25 and 0.7 GeV/c^. The Kolmogorov-Smirnov 
test gives a maximum deviation of 0.1% which is small 
compared to the uncertainty of the absolute normaliza- 
tion of the mixed-event background. The corresponding 
Kolmogorov-Smirnov p-value is also shown in the right 
panel of Fig. |l^ as a function of mass. It is small below 
0.55 GeV/c^, where there is some contribution of the cor- 
related background. However, for rriiow > 0.55 GeV/c^, 
the p-value is ^90%, therefore confirming that the hy- 
pothesis of compatibility of the two distributions for 
TOiow > 0.55 GeV/c^ is valid for any commonly used sig- 
nificance level. 

Figure ^ shows the like-sign mass distribution for 
real and mixed events in the different centrality bins 
used in the analysis. The bottom panels show the ra- 
tio (iV±± - B§°|^'=)/B5°|^'^ which are fit to a constant for 
mee > 0.55 GeV/c^. The fit results for all ccntralities are 
reported in Table III together with the results of the sta- 
tistical tests described above. The results reported in Ta- 
ble [II demonstrate that the agreement between real and 
mixed-event like-sign mass spectra demonstrated for min- 
imum bias collisions also holds for all centrality classes. 

Figure |l^ shows the like-sign mass distribution for real 
and mixed events in different pt bins. The ratio (A^±± — 
^comby^comb gj^Q^^^g g, good agreement between real and 
mixed-event distribution for all bins. This demonstrates 
that, within the statistical error of the foreground, there 
is no deviation from uncorrelated combinatorial behavior 
for masses above 0.55 GeV/c^ in any range. 

The absolute normalization of the unlike-sign combi- 
natorial background is determined with an uncertainty 
given by the statistical error of the measured like-sign 



spectra for rriee > 0.7 GeV/c^ of 0.12%. The mass inter- 
val (mee > 0.7 GeV/c^) is sufficiently large to achieve the 
desired statistical accuracy and is conservatively chosen 
to exclude any possible region which may be contam- 
inated by the correlated background. The results are 
stable when varying the lower end of the normalization 
region between 0.55 and 0.7 GeV/c^, and the difference 
is included in the systematic uncertainty. 

Because the pair cuts re move m ore like-sign than 
unlike-sign pairs, (see Section HID 1 ) the normalization 
factor is corrected by this asymmetry, which is estimated 
with mixed events to be 1.004 ± 0.002, independent of 
centrality. 

For the various centrality bins and minimum bias col- 
lisions, the 0.2% uncertainty on the event rejection is 
added in quadrature to the uncertainty on the normal- 
ization, which is determined by the statistics of the like- 
sign pairs. Since the ratios shown in Figs. ^, |l^, and 
show no systematic deviation in shape between the 
like-sign real- and mixed-event distribution for rriee > 
0.55 GeV/c^, the uncertainty due to the shape is negli- 
gible compared to the uncertainty on the normalization. 
The total uncertainty on the combinatorial background 



is given in Table IV 



TABLE IV: Systematic uncertainty on the combinatorial 
background for p+p and various centrality classes of Au ± Au 
collisions. 



Centrality 



^gcomb^^comb ^^^^ 



0-92% 
0-10% 
10-20% 
20-40% 
40-60% 
60-92% 
p±p 



0.25 
0.25 
0.30 
0.35 
0.85 
2.50 
3.00 
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FIG. 11: (Color online) Like-sign distribution for real N±± 
and mixed events B™™'' for different centrality data sets. The 
bottom panels show the ratios (iV±± - B'l^^^) / B'^^'° with 
different scales. The left ones show all the mass range and 
allow to identify the correlated background at low masses. 
The right ones focus on mass region where we fit. 



FIG. 12: (Color online) Like-sign distribution for real and 
mixed events for different pr bins. The bottom panels show 
the ratio of (7V±± - B'^^'°) / B'^^^ with different scales. The 
left ones show all the mass range and allow to identify the 
correlated background at low masses. The right ones focus 
on mass region where we fit. 



This translates into a systematic uncertainty SS on 
the signal S of SS/S = (5B™"*/B'=°'"'')/(S'/B™'"'^). 
jljcomb^^comb j^^j. jnijijjn^jn j^jag collisions and all cen- 
trality bins are listed in Table IV. Figure |l^ shows 
gjj^corah different ranges of pair-p^ and centrality. 



Correlated Background After subtracting the com- 
binatorial background, the remaining pair distributions, 
like- and unlike-sign, are considered correlated pairs, 
where the like-sign distribution only contains correlated 
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FIG. 13: (Color online) Ratio of signal to background for different ranges of pair-pr (left) and centrality classes (right). 



background pairs while the unhke-sign also contains the 
signal. 

N±± - B^^^^ = (13) 
N+- - B;™* = + S+- (14) 

The correlated background B'^°" arises from two sources. 

The first source is "cross pairs". They occur if 
there are two e+e~ pairs in the final state of a meson, 
e.g. double Dalitz decays, Dalitz decays followed by a 
conversion of the decay photon or two-photon decays 
followed by conversion of both photons. Besides the real 
unlike-sign signal, this leads to like- and unlike-sign cross 
pairs. While all mesons in principle produce cross pairs, 
all contributions above the r] mass (0.55 GeV/c^) can be 
safely neglected. Like- and unlike-sign cross pairs were 
simulated using our hadron decay generator including 
the PHENIX acceptance (Equation ^). Because their 
rate is proportional to the Dalitz decay probability 
only, the ratio of cross pairs to pions is independent of 
centrality. 

The second source is "jet pairs" . They are produced by 
two independent hadron decays yielding electron pairs, 
either within the same jet or in back-to-back jets. Jet 
pairs were simulated using the minimum bias events gen- 
erated with PYTHIA as described above. As noted 
previously, correlated pairs from the same jet typically 
have small mass and large pr while those from back-to- 
back jets have large mass and smaller pr- 

Correlated background pairs equally populate like- 
and unlike-sign combinations. Since the like-sign spec- 
trum measures only the background, we can determine 
the cross and jet pair yields by simultaneously fitting 
simulated cross and jet pair distributions to the mea- 
sured correlated like-sign pairs, after subtraction of 
combinatorial background. The resulting normalization 
factors, one for cross-, one for jet-pairs, are then applied 
to the unlike-sign correlated background. 



Correlated Background in p + p data 
Panel a) of Fig. |l^ shows like-sign pair distributions: 
^^±±, 5±±^ and B^^ = N±± - B±±. Panel b) shows 

analogous distributions for unlike-sign pairs: 7V_| , B^ 

and B^j^!_' = 7V+_ - B+_. 

The like-sign correlated background distribution is fit 

to 

B™J5 = A ■ cross + B ■ jet (15) 

The resulting normalization factors, A and B, are then 
applied to the unlike-sign correlated background. The 
unlike-sign signal results from subtracting the correlated 
background from all correlated unlike-sign pairs. 

Correlated Background in Au + Au data 

In Au + Au data the like-sign mixed-event distribution 
reproduces the mass dependence of the real-event dis- 
tribution not only in region of Equation but for all 
masses above 0.55 GeV/c^ for every centrality, as shown 
in Fig. and at every px, as shown in Fig. |lj. This 
means that there is no room for correlated background 
for ruee > 0.55 GeV/c^. Here the contribution that typ- 
ically arises from back-to-back jets in the Au -I- Au data 
is indeed expected to be different than p + p because of 
the observed jet modifications pOf . 

We therefore separate the jet distribution into "near- 
side" {jetnear- < tt/2) and "away-side" contributions 
{jetaway- > 7r/2) and we fit the hke-sign correlated 
background distribution to the sum of 

B'^" = A ■ cross + B ■ jetnear + C ■ jetaway (16) 

Panels c) and d) of Fig. |lj show the like- and unlike- 
sign pair distributions, the normalized mixed-event back- 
ground and the distributions after subtraction. We note 
that the like-sign yield is well described by the sum of 
combinatorial and correlated background, and that the 
contribution from "away-side-jet pairs" is consistent with 
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FIG. 14: (Color online) Raw dielectron spectra in p+p (left) and An + Au (right) collisions. The top panels show like-sign pairs 
N±± as measured in the experiment, the combinatorial background from mixed-events B^j?™'', the correlated pair background 
B±± obtained by subtracting the combinatorial background, and the individual contributions from cross and jet pairs to the 
correlated background (see text) . The bottom panels show the same distributions for unlike-sign pairs. The correlated like-sign 
background B±± , is normalized to the measured like-sign pairs remaining after subtracting the combinatorial background 
N±± — B™™'', and the same factors are applied to the unlike-sign distribution B™"'. 



zero, i.e. C = 0, as the mixed-event distribution was nor- 
malized to the real data in the IMR: the "away-side-jet 
pairs" are therefore not shown. The unlike-sign signal 

(S^ ) is obtained by subtracting from the distribution 

of all pairs the mixed-event combinatorial background 
^^comb-j correlated background {B^") normal- 

ized with the factors A,B and C measured in the like-sign 
spectrum. 



3. Like-sign Subtraction Method 

The event-mixing technique, along with the determi- 
nation of correlated backgrounds, used to estimate the 
background contribution to the measured dilepton mass 
spectra, was developed in order to get around the prob- 
lems introduced in a traditional like-sign background cal- 
culation by the asymmetric PHENIX acceptance. With a 
like-sign calculation we need to make only the assumption 
that the mass dependence of the correlated background is 
symmetric for like- and unlike-sign pairs. We do not need 
to make any assumption about the decomposition of the 
background into "cross", "jet-near-side" and "jet-away- 
side" pairs. The like-sign distribution measures only the 
background. Correlated and uncorrelated backgrounds 
are charge symmetric, i.e. they yield the same number 
of like-and unlike-sign pairs, with the same distribution. 



The measured distributions however are different because 
of the acceptance. 

The acceptance difference between like- and 
unlike-sign pairs can be measured with the ratio 

of mixed-event distributions, i.e. 



where B^°^^ , B^^"*^ , B^i"'^ 



2^Bco„,b.Bco.„b 

are the number of 
e ' e , e ' e ' , e e pairs in a given mass and pt bin. 
The like-sign distribution can therefore be acceptance- 
corrected and then subtracted from the unlike-sign 
pairs. 



B 



comb 



2^^^comb _ ^comb 



(17) 



This technique however measures the background distri- 
bution with similar statistical error as the foreground, 
therefore the resulting error on the signal is much larger 
than when using the mixed events technique. 

One can also use a "hybrid" method where first the 
mixed events are normalized in the region where real and 
mixed-event distributions agree (region in Equation |l^ 
ioY p + p and lUee > 0.7 GeV/c^ for Au -|- Au). Then the 
correlated background is obtained by correcting the sub- 
tracted like-sign distribution by the acceptance difference 
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FIG. 15: (Color online) Unlike-sign dielectron spectra obtained by subtracting the background (combinatorial and correlated) 
with the methods explained in the text for p + p (left) and Au + Au (right). The RMS/MEAN of all the spectra shown in the 
bottom panels allows to assign a systematic uncertainty on the signal due to correlated background subtraction. 



4- Systematic Uncertainty of the Background Subtraction 



> 0.55 GeV/c^; 



The background subtraction is the most critical aspect 
of this analysis, therefore it is crucial to assign a proper 
systematic u ncertain ty to it. 

In Section [II D 2| we reported the systematic uncer- 
tainty on the normalization of the combinatorial back- 
ground determined by the statistical accuracy of the 
measured like-sign yield (see Table IV ) . The uncertainty 
on the shape of combinatorial background is everywhere 
negligible compared to the uncertainty on the normaliza- 
tion. 

To evaluate the systematic uncertainty due to the cor- 
related pairs we compare the results obtained with the 
two subtraction methods. It is important to note that the 
contribution of the correlated background is everywhere 
small compared to the signal. Therefore any uncertainty 
in its estimate would result in a small uncertainty on the 
signal. 

In the first subtraction method the correlated back- 
ground B'i^" was calculated using Monte Carlo simula- 
tions. We have estimated potential spectral shape mod- 
ifications of B"^" in several ways, shown in Fig. |l^. In 
Au + Au 

• we have not subtracted any jet background (this 
would clearly result in an upper limit); 



• we have subtracted both near and away-side using 
for the away-side the same normalization as for the 
near-side (this would clearly result in a lower limit); 

• we have subtracted the near side normalizing it for 
70° <A(t)< 110°; 

• we have subtracted the near side normalizing it for 
70° < A(j) < 110° and we have corrected the shape 
for the elD efficiency. 

In the second subtraction method, both in p + p and 
Au -f- Au we use the like-sign corrected for the different 
acceptance, or a hybrid method which uses the mixed- 
event distribution for the combinatorial background and 
the like-sign distribution corrected for the different ac- 
ceptance for the correlated background only. 

Since the acceptance for pairs is a function of mass and 
Pt, we have checked that for different e+e^ pair sources, 
which span reasonable variations in mass and px shapes 
of the e+e^ pairs, the relative acceptance is unchanged. 
For this purpose in p + p we have calculated the relative 
acceptance: 

• using data from p + p run only; 

• using data from Au + Au run only; 



• we have subtracted the near side normalizing it for 



• using data from pythia real events; 
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• using data from PYTHIA mixed events. 

We added the results from these two methods in Fig. |l^ 
as well. The agreement within the different methods (see 
RMS/MEAN in the bottom panel) allows to assign a 
systematic uncertainty on the signal. 



E. Raw Mass Spectrum 



+e pairs, 



Figure |l^ shows the mass distribution of e 
the normalized mixed-event background {B), and the sig- 
nal yield (S) obtained by subtracting the mixed-event 
background and the correlated background (cross and 
jet pairs) for Min. Bias Au + Au collisions. The right 
panel shows the signal-to-background ratio {S/ B'^°"'^^). 
The systematic uncertainties (boxes) reflect the uncer- 
tainty from the background subtraction, which is given 
by Ss/S = 0.25% • B'=°"'^/S, added in quadrature to 
the uncertainty due to the correlated background sub- 
traction, conservatively assumed to be around 10%S be- 
low 0.55 GeV/c^. Despite the smaU S'/B™"'^ ratio, an 
e+e" pair continuum is visible up to 4.5 GeV/c^. Due to 
the limited statistical precision in Au + Au binning was 
chosen such that the bins near the (f> and uj meson are 
centered at the nominal meson mass. The bin width was 
chosen to correspond approximately to twice the mass 
resolution observed in p + p collisions. In this way the 
and Lo mesons can be resolved from the continuum. 



F. Runs with Increased Conversion Material 

In order to check the background subtraction, a sub- 
set of data (5x10^ events) was collected with a brass 
sheet of 1.68% radiation length {Xq) wrapped around the 
beam pipe to increase the number of photon conversions. 
We make an estimate on the uncertainty in the radiation 
length (5%) by comparing photon conversions in data 
and simulation Because the additional conversion 
leads to an increased electron multiplicity, in this data set 
the combinatorial and correlated background (_B) contri- 
bution is larger by a factor of ^2.5. 

If there is a systematic bias in background normaliza- 
tion, and the yield of background is off by a small frac- 
tion /, it can lead to significant difference between the 
"apparent" signal and the "true" signal S'*'""'^. The dif- 
ference should be larger in the converter run since signal- 
to-background is smaller. Thus the "apparent" signal 
in the converter run should become larger than the "ap- 
parent" signal Snc in normal run without the converter. 
The relation between 5*™°, Sc, Sj^c can be written as: 



5nc = iVNC - (1 - /) • S 

Sc ^ Nc- 2.5 • (1 - /) • B 
= 5'*™'= -I- 2.5/B 



(18) 
(19) 



where Nq (-/Vnc) is the foreground of all e+e" pairs in 
converter (non-converter) runs. If we divide Equation |l^ 
by B, we obtain: 



NC 



IB 



B B B 

gtruc ^ JB= iS/B)B 



(20) 



Then the ratio between the apparent signal in converter 
and non-converter runs is given by: 



5, 



^.truc _^ 2.5fB 



NC 



^true 
^.truc ^ 



fB + 1.5fB 



{S/B)B + 1.5,fB 

(S/BjB 
S/B + 1.5f 



S/B 



Therefore 



/ 



'NC 



B L5 



(21) 



(22) 



This means we can use the agreement between con- 
verter and non-converter runs 5'c/>S'nci shown in Fig. [l^ , 
and the signal-to-background ratio S/ B also shown in the 
insert of Fig. |l6|, to constrain a potential bias / in the 
background normalization. 

In 0.3 < mee < 0.75 MeV/c^ 



Sc_ 

Snc 



1.29 ±0.92 



and 



— = (7.47 ±0.55)10-3 
B 

constrain / = (0.14 ± 0.46)10~^. If we extend the mass 
range to lower values 0.04 < rriee < 0.75 MeV/c^, 



Sc_ 

Snc 



1.08 ±0.28 



and 



4 = (2.75±0.05)10-^ 
B 

the constraint becomes / = (0.15 ± 0.51)10^^. 

The agreement between the converter run and the nor- 
mal run confirms that the systematic bias in the back- 
ground normalization is small (0.15±0.51%). The result 
is consistent with our estimate of the systematic un- 
certainty in the background normalization (0.25%). It 
should be noted that the converter run provides an inde- 
pendent test on the background normalization. 
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FIG. 16: (Color online) Uncorrected mass spectra of all e'^e~ pairs, mixed-event background (5"=°™'') and signal (S) in 
minimum bias Au + Au collisions. Statistical (bars) and systematic (boxes) uncertainties are shown separately. The signal 
from the runs with additional converter material is shown with statistical errors only. The two spectra are normalized by the 
number of events. The insert shows the S/B""'"^^ ratio. The mass range covered by each data point is given by horizontal bars. 



G. EfHciency Correction 



of the pair cuts (epair and ep^"''*) is taken into account 



pair 



The e^e" mass spectra are corrected for the total pair 
reconstruction efficiency e^°^^^ to give the e+e~ pair yield 
inside the PHENIX aperture as defined by Equation p|: 



total cID _ live 

pair pair pair 



,0V 



ghost 



pair pair 



^pair 



(for Au + Au) 
(for p + p) 



ERT 
^pair 



(25) 



dN 
dmp„ 



1 



Ncvt Amee e 



total 
pair 



(23) 



The pt spectra are further corrected for the pair geo- 



metric acceptance (cp^j.) 
one unit of rapidity: 



to give the e^e -pair yield in 



1 



1 



total ,goo 
pair 



(24) 



The p + p data are further corrected by the factor 
l^fi^ where ebbc = 54.5 ± 5% is the BBC efficiency and 
ibias = 7 9 lb 2% is the BBC trigger bias, described above 
(Section pIAl). 



The total pair reconstruction efficiency 6*°'^ depends 
on the single electron efficiency for reconstruction and 
elD (Esi™!;,), the efficiency from the detector live area 

^livc 



V single 



the occupancy efficiency (e°in'gic) (for Au -I- Au). 
The occupancy efficiency (e°ingio) is different for each 
Au -I- Au centrality class and is described in Section III J| . 
The p + p data are also corrected for the ERT trigger ef- 
ficiency, described in Section III I. In addition, the effect 



These efficiencies depend on the elD cuts used to de- 
termine the electron sample, therefore they do factorize 
only on the condition that the electron sample used to 
calculate them is the same for all them. 

The pair-elD and reconstruction efficiency 6°^^? and 

pair-live efficiency ep™r' well as the ERT efficiency 

ERT 
pair 



foT p + p data are derived, as a function of mass 
and pair-pT, from the corresponding single electron ef- 
ficiency (eJi^gie, e^^gie and e^JJ^^ respectively) using the 
pair kinematic properties implemented in our hadron de- 
cay generator, as explained below. 

is the fraction of signal loss due to track 



The 



oID 
single 



reconstruction and elD cuts within the detector active 
area. It depends only on the momentum of the track. 
The shape is very similar in p + p and Au -I- Au, but the 
scale is different since different cut values are used. It is 
calculated as follows: 



single (Pt) 



llvf/dp^' 



(26) 



The dN™/dp?p is the p?p distribution of the input elec- 
tron yield that falls into the real PHENIX acceptance, 
which includes the boundary described by Equation ^ 
and the active areas of the detector. dN°'^^ / dpfp is the 
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distribution of the output electron yield in the same 
acceptance after passing all the elD cuts. 

The input distribution comes from a simulation of 
450M 7r"'s fiat in phase space {0 < pr < 25 GeV/c, 
\y\ < 0.5, and < (j) < 27r) with the branching ratio of 
the 77° Dalitz decay set to 100% to enhance the sample 
of e~^e~ pairs per event. These events were processed 
by the full GEANT simulation program of the PHENIX 
detector [|2) that includes the details of the detector re- 
sponse. The output simulation data files were processed 
by the event reconstruction chain of PHENIX. Standard 
elD cuts are apphed to the output. The reconstructed 
pfn of each output electron is weighted according to 
with a realistic exponential px weight. 

Figure shows the single electron (positron) efficiency 
as a function of p^ for p + p and Au + Au. The different 
scales corresponding to the two data sets are due to more 
stringent cuts applied in the Au + Au data set. The 
band around the curve shows only the p^ dependence 
of the systematic uncertainty corresponding to a shift 
of ± 0.1 GeV/c of the efficiency curve. This in turns 
leads to a distortion of the e+e" mass shape shown by 
the band in Fig. |l^. The total uncertainties on the pair 
reconstruction, including the range of applicability, are 
reported in Table VI. 

The Cs'ingie is the fraction of signal loss due to inactive 
areas of the detector. The active areas of the detector 
are parameterized as a function of the particle momenta 
and azimuth using real data. There are small differences 
between p + p and Au + Au. 

In the EXODUS cocktail of hadron decays we have im- 
plemented the parameterization of the single electron 
(positron) efficiency (shown as a curve in Fig. |l^. The 
single electron (positron) efficiency was applied as a 
weight to each track. The pair will therefore get a weight 
given by the product of the electron and the positron 
weight. This weight, a function of mass and pair-p^, 



represents e 



cID 



To calculate e 



live 



a fiducial cut corre- 



^pair" ^^^^^^^^^ ^pair' 

sponding to the detector active areas is also implemented 
in EXODUS and a pair is rejected if at least one track falls 
out of the active areas. We therefore determine the prod- 



uct e! 



!ID 



X 



double differentially as a function of mass 



X el^=fr as a function of mass. 



and pair-pT 

Figure |l^ shows e' 
At high masses the efficiency is constant with a value 
nearly equal to the square of the single electron efficiency. 
At low mass the pair efficiency results from the con- 
volution of the single electron reconstruction efficiency, 
which drops toward low momentum, and the geometric 
acceptance, which effectively truncates the single elec- 
tron pt distribution {p^ > 0.2 GeV/c). For pairs with 
0.4 < niee < 0.8 GeV/c^ the efficiency drops as a conse- 
quence of the drop at \ow-pT of the single electron effi- 
ciency. However, for niee < 0.4 GeV/c^, the lower limit 
on single electron px results in a larger average momen- 
tum; and the pair-efficiency consequently increases. The 
band around the 6°™. x e^^f^ curve in Fig. |l^ again shows 
only the possible distortion of the mass distribution, due 



to the uncertainty shown in Fig. |T^. It reaches a maxi- 
mum of 8% around mee = 0.4 GeV/c^. 

The £^1°^^ represents the loss of real pairs which acci- 
dentally fulfill the detector overlap criteria. It is deter- 
mined by the corresponding loss in mixed events as: 



ghost ^ dBll^/jdrrieedpT) 
'P^"- dB2lJ{dmeedpT) 



(27) 



where dBl {dB2) are the mixed-event unlike-sign pair 
(mass, px) distribution with (without) applying the over- 
lap pair cuts. 

The CpJ;^ represents the loss of real pairs which are 
accidentally oriented like conversion pairs in the mag- 
netic field. At low masses it is calculated using the tt*^ 
GEANT simulations described above. Above the pion 
mass, we use our hadron decay generator exodus where 
we have implemented an empirical smearing for the de- 
tector resolution in the determination of the magnitude 
and direction of the momentum vector. 

Figure |l^ shows the total pair-efficiency as a func- 
tion of invariant mass for different ranges of pair-p^ for 
Au + Au collisions. 

The systematic uncertainty on e'^^^ x e^^lf^. is given by 
twice the uncertainty on the single-particle efficiency, de- 
termined by varying the cut values. In addition, un- 
certainties in the active areas have been determined by 
varying the active areas in the simulations. 

The systematic uncertainties depend somewhat on pT 
and mass, they are largest at low pr and low mass where 
they reach 14.4% for p+p and 13.4% for Au -|- Au. We use 
these errors as an estimate for all momenta and masses. 
Since the pair-efficiency was derived using realistic kine- 
matics, the uncertainty arises mostly from low-p^ tracks 
and neglecting the mass and pair-py dependence is a con- 
servative approach. 



The 



ghost 



and et^. affect only the LMR. Therefore, 

pcill [J (±11 ' 

the corresponding systematic uncertainties have been in- 
dependently evaluated by varying the cut values. We 
estimate an uncertainty of 5% due to the ghost cut and 
6% due to the photon conversion cut on the final e+e^ 
yield and we apply it to the LMR by varying the cut 
value in a suitable range. 

The (mass, pt) distributions from every individual 
centrality bin has been corrected using the same 2D ef- 
ficiency correction function calculated for Min. Bias. 
This procedure avoids an assumption for the input kine- 
matic distribution. However, since the individual cen- 
trality bins have limited statistics compared to the Min. 
Bias sample, this procedure may suffer statistical fluctu- 
ations. Alternatively we use an effective efficiency correc- 
tion, as a function of mass only, obtained by weighting 
the 2-dimensional (mass, px) corrections with a realistic 
Pt distribution provided by the minimum bias data set. 
This curve is given by the product of all the curves shown 
in Fig. |l8|. We take the difference of 10% between the ID 
and the 2D corrected result as an additional systematic 
uncertainty for the efficiency correction of the different 
centrality data sets (see Fig. p9|). 
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FIG. 17: (Color online) Single electron (positron) efficiency as a function of for electrons (left) and positrons (right) in 
p + p (top) and An + Au (bottom) collisions. 
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FIG. 18: (Color online) Different components of the total pair efficiency as a function of pair mass, obtained with the procedure 
described in the text, for the p + p and the Au + Au collisions. 



We verified this efficiency using a GEANT simulation 
of 1 million e+e- pairs. Roughly half were generated flat 
in mass (0-4 GeV/c^), (0-4 GeV/c), azimuthal angle 
(0-27r), and rapidity {\y\ < 0.5). The other half were gen- 
erated with a probability inversely proportional to px in 



order to enhance the statistics in the low-mass and low- 
Pt region, where the efficiency varies most. Only pairs 
with both an electron and a positron in the ideal accep- 
tance (given by Equation H) are processed by GEANT 
and reconstructed with the same analysis chain. The 
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FIG. 19: (Color online) ep^^^ x e^^^ as a function of invariant 



efficiency is determined double differentially in px and 
mass of the e+e~ pair. This second method gives consis- 
tent results, but is limited by MC statistics. It was used 
as cross-check in the final analysis. 



0.1 GeV/c^ wide mass bins. The cocktail consists of the 
sum of polarized Dalitz decays and unpolarized vector 
mesons, while the flat simulations are always unpolar- 
ized. The shape of the acceptance is very similar, and 
the relative normalization agrees within 5% in the lowest 
mass bins, and better for higher bins. Based upon this 
comparison, illustrated by the ratio in the bottom figure, 
we assign an upper limit of 10% (marked by the lines in 
the figure) for the systematic uncertainty of e^°^j.. 

Since it arises from the independent fragmentation of 
two charm quarks, the contribution of charmed meson 
decays has a different acceptance. This component has 
been simulated with PYTHIA | 77 and normalized accord- 
ing to the cross section measured in [|6| scaled by A'^coU- 
However, due to the observed modifications of charm 
quarks in the medium j^], the acceptance could poten- 
tially be different than what is simulated by pythia. For 
niee < 0.5 GeV/c^ the charm contribution is negligible. 
For 0.5 < TOee < 1 GeV/c^ a systematic uncertainty of 
5% due to the uncertainty of the charm cross section 
(cr^g = Nc^ii X 567 ± 57"*=^* ± 224">"'' /ib §, |6)) has been 
added in quadrature to the other systematic uncertain- 
ties on efZZ- 



H. Acceptance Correction 



I. Trigger Efficiency (p -I- p) 



In addition to the efficiency corrections, the pr spectra 
are also corrected for the detector geometric acceptance 
^palr' e"*'e~-pair yield over the full azimuth 

in one unit of rapidity. ep°°i. accounts for the fraction 
of pairs produced in one unit of rapidity over the full 
azimuthal range that are lost because either one or both 
particles from the pair miss the PHENIX detector 



- (in PHENIX acceptance) 



pair 



dp^drue 



:(|yee| <0.5) 



(28) 



pair 



has been calculated using a Monte Carlo simula- 
tion of Dalitz decays of pseudoscalar mesons (tt*^, yy, 77') 
and direct decays of vector mesons {p,uj,(j), 3/ip, 
For all mesons the rapidity distribution is assumed to 
be flat around mid-rapidity. This assumption is well jus- 
tified as PHENIX measures in \yee\ < 0.35, where the 
natural distribution is flat. The acceptance is therefore 
uniform in \yee\ < 0.35 and we do not assign a system- 
atic uncertainty to it. The px distributions are taken 
from PHENIX measurements (see Section [V), and me- 
son polarizations are taken from p3| . The acceptance 
correction is performed double differentially in pair-mass 
and pt with 0.005 GeV/c^ bins in mass and 0.1 GeV/c 
bins in px- 

The systematic uncertainty due to the px parameteri- 
zation and the polarization of the Dalitz pairs is studied 
with a simulation of unpolarized pairs with a flat dis- 
tribution in mass and pair-py- In Fig- ^ we compare 
the acceptance for the full hadronic cocktail and the sim- 
ulation of unpolarized pairs as a function of pair-p^ in 



The efficiency of the ERT trigger (e™^) 'm p + p col- 
lisions, as a function of pair mass and px, is determined 
with a fast Monte Carlo simulation of pairs which uses 
a parameterization of the single electron ERT efficiency 
^ingic' "^^^ single electron ERT efficiency is determined 
using the Min. Bias data set. The online level-1 trig- 
ger decision is recorded in the Min. Bias data set even 
though it is not used to select the event. We require that 
the trigger tile that fires the ERT trigger is used by the 
electron candidate selected by the offline analysis. The 
ratio of triggered electrons relative to all electrons candi- 
dates gives the trigger efficiency as function of the single 
electron and is shown in Fig. The trigger efficiency 
at the plateau is « 60%, consistent with the fraction of 
the active trigger tiles. The trigger part and the offline 
part of the RICH and EMCAL read-out are handled in 
separate electronics chains, and the trigger part has more 
noisy tiles which are masked out. This results in less ac- 
tive areas in the trigger. The efficiency is fit to the sum 
of two Fermi functions 



fiPT) 



£0 ■ 9{px - 0.5) 

g-(pT-po)/fc -|- 1 



3 ■ 9{px + 0.5) 

-iPT-p'o)/k' ^ I 



(29) 



where ep (e'o), Po (p'o) and k {k') are free fit parameters 
and 6 is the usual Heaviside theta function. 

This parameterization is used in a fast Monte Carlo 
simulation of hadron decays into e+e~. Each electron 
fires the trigger with a probability given by the trigger 
efficiency. We require that both electrons are within the 
PHENIX acceptance and at least one of them fires the 
trigger. The mass and px distribution of those pairs is 
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FIG. 20: (Color online) e^e~-pair acceptance as a function of pair-pr for different mass ranges for the hadronic cocktail (line) 
and a flat distribution of e"'"e~ pairs (empty circles) and ratio of the two (bottom figure). The lines mark the limits of the 
systematic uncertainty. 



compared to the one obtained with no trigger require- 
ment. The ratio is the pair trigger efficiency, and is shown 
in Fig. ^ as a function of e+e^ invariant mass. The 
structures in Fig. ^ result from the turn-on curve of the 
trigger threshold convoluted with the acceptance of the 
detector. 



The systematic uncertainty has been studied by vary- 
ing the parameterization of the single electron ERT effi- 
ciency, as well as varying the elD cuts which define the 
reference sample. The effect of changing the elD cuts 
leads to a larger uncertainty at low-pr and low-masses. 
We therefore assign a systematic uncertainty of 20% for 
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FIG. 21: (Color online) Trigger efficiency for single electrons 
Single ^ function of pr of the ERT trigger in p + p collisions 
determined from the Min. Bias data set. 
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FIG. 22: (Color online) Trigger efficiency e^Jf^ for e+e -pairs 
as a function of e^e~-pair invariant mass. 



rriT < 1 GeV/c^ and 5% elsewhere. Finally, the system- 
atic uncertainty includes potential shape distortion to the 
efficiency correction due to the variation of active trigger 
tiles during the data collection. 

Fig. ^ compares the invariant mass spectra for the 
p+p data obtained with the Min. Bias and the ERT 
data sets. The ERT data set has been corrected by ERT 
trigger efficiency e^^J and the total pair reconstruction 
efficiency ep°*^'i while the Min. Bias is corrected only 
for the 6*°*^'- This comparison of the two data sets con- 
firms that the Min. Bias and the ERT agree within their 
respective statistical errors. 
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FIG. 23: (Color online) Invariant mass spectra for p + p data 
with the Min. Bias (hollow) and ERT (solid) data sets. The 
agreement between the two data sets is excellent. 



detector hits from real Au + Au events. Next, these new 
files containing the embedded are run through the 
entire reconstruction software. As the particle density 
reduces the efficiency but does not introduce additional 
Pt dependence, the occupancy correction can be factored 
out. Since all the detectors used in the analysis are lo- 
cated after the pair has been opened by the magnetic 
field, the pair embedding efficiency in each centrality bin 
is defined as the square of the single electron embedding 
efficiency 



/ occ \2 
V^singlc/ 

# reconstructed from embedded data 
# reconstructed from single track data 



where a reconstructed particle from embedded data is 
required to have most of its detector hits associated with 
hits from the simulated particle. 



TABLE V: Embedding efficiency for different centrality 
classes of Au + Au collisions. 



Centrality 



^sing;lc 



0-10% 
10-20% 
20-40% 
40-60% 
60-92% 



0.86 
0.91 
0.93 
0.97 

0.99 



0.74 
0.83 
0.87 
0.95 
0.98 



J. Occupancy Correction (Au + Au) 

In Au + Au collisions there is an additional efficiency 
loss of particle detection due to the presence of other 
particles nearby. To study this effect single electrons and 
positrons are simulated through the GEANT simulator 
of PHENIX and then embedded into data files containing 



Table |V| displays the embedding efficiencies for the cen- 
trality classes used in the analysis. For the minimum bias 
results we have weighted the occupancy correction by the 
fraction of pairs in each centrality class. Since most of 
the yield is concentrated in the most central classes, the 
resulting pair efficiency is 0.81 instead of the square of 
the minimum bias value 0.96^ = 0.92. 
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A second, data-driven method was developed to de- 
termine 6°!^;^.. This method uses the conversion pairs cre- 
ated at the beam-pipe to select a track sample of pure e+ 
or e~ from real data. The conversion pairs originating 
at the beam-pipe result in a clear invariant mass peak 
around ^ 20 MeV/c^ (see Fig. ^). We assign tight 
elD cuts to only one track of a pair and measure the ef- 
ficiency of the other track of the same pair without any 
elD cut applied. Differences between the two methods 
are accounted for in a 3% systematic uncertainty pl|. 



Most of these uncertainties are mass-pr correlated, i.e. 
all points move in the same direction, but not by the same 
factor. Only the BBC and trigger bias (in p + p) and the 
occupancy (in Au + Au) are normalization uncertainties 
in which all points move by the same factor independent 
of mass and px- Since those uncertainties are small com- 
pared to the total uncertainty, they are included in the 
total uncertainty, without plotting it separately. 

IV. COCKTAIL OF HADRONIC SOURCES 



K. Systematic Uncertainty 

The systematic uncertainties are summarized in Ta- 
ble VI. The uncertainties in Table VI are categorized by 
type: 

• A: point-to-point uncertainty uncorrelated between 
mass or px bins, 

• B: correlated uncertainty, all points move in the 
same direction but not by the same factor, 

• C: an overall normalization uncertainty in which 
all points move by the same factor independent of 
mass or px- 

The uncertainty on pair reconstruction efhciency includes 
elD cuts, geometric acceptance, and run-by-run fluctu- 
ations. The uncertainty due to the conversion rejection 
cut and the overlap cuts as well as the uncertainty due to 
the ERT and minimum bias trigger efficiency (for p + p) 
and occupancy (for Au -I- Au) are also listed. These 
uncertainties are included in the final systematic uncer- 
tainty on the invariant e"'"e~-pair yield. The uncertain- 
ties deriving from reconstruction and occupancy do not 
have a strong px dependence, so we keep the assigned 
values for every px bin. The pair cuts and the conver- 
sion rejection are localized in mass (m < 0.6 GeV/c^) 
and are rather py-independent. The uncertainty on the 
combinatorial background is the largest contribution to 
the systematic uncertainty in Au -t- Au and the value is 
estimated to be 0.25% -B/S. S/B rises with px, there- 
fore the uncertainty on the combinatorial background has 
been propagated separately for each px bin. The uncer- 
tainty on the correlated background is approximately 2- 
3% in p + p and 10% in Au + Au for uiee < 0.6 GeV/c^ 
and pT-independent. However it increases in p + p to- 
wards high masses. In the p + p data the uncertainty 
on the ERT trigger efficiency, as well as on the BBC and 
trigger bias is included. For individual centrality bins, we 
add 10% uncertainty arising from the px dependence of 
the efficiency correction: this was obtained from the dif- 
ference between the ID-corrected mass spectra and the 
2D-corrected mass spectra. For the px spectra, which are 
further corrected by the geometric acceptance, we added 
10% uncertainty from the acceptance correction and 5% 
for ruee > 0.5 GeV/c^ due to the charm contribution. 



In this Section we describe the methods used to calcu- 
late the pair yield expected from hadronic decays which 
will be compared to the experimental data. 

We model the e^e~ pair contributions from hadron 
decays using the decay generator exodus, exodus is a 
phenomenological event generator that allows to simulate 
the phase-space distribution of all relevant sources of elec- 
trons and electron pairs and the decay of these sources. 
Also it allows to include the filtering for the geometrical 
acceptance and the detector resolution. The relevant pri- 
mary mesons that involve electrons in the final state are 
tt'^, rj, rj' , p, uj, (j>, J/tjj and ip' . 

We assume that all hadrons have a constant rapidity 
density in the range |A7/| < 0.35 and a uniform distri- 
bution in azimuthal angle. Transverse momentum distri- 
butions are largely based on measurements in PHENIX. 
The key input is the rapidity density dN/dy of neutral 
pions, which we determine from a fit to PHENIX data on 
charged and neutral pions in H, |6^, |l| with a modified 
Hagedorn function, given by: 



E 



dp'^ 



(30) 



Fit parameters and dN/ dy for p + p and Au + Au are 
given in Table VII. 

For all other mesons we assume mx scaling, replac- 
ing Px by m? — rn^ + (p^/c)^, where m is the mass of 
the meson and we fit a normalization factor to PHENIX 
data g, H, g, I^H, H, H, ^, ^, H, where 
available. Figure p4| shows the excellent agreement with 
published PHENIX data. The rj meson is measured only 
at higher px] however, in the p + p collisions the fit is in 
good agreement with the px distribution of kaons, which 
have similar mass (see discussion below). 

In order to extract the meson yield we integrate the 
fits over all px- Results, systematic uncer tainties, and 
references to data are given in Table VIII and the ra- 
tio of the integrated yields meson /i:^ are compared for 
p + p and Au + Au data. For the p meson we assume 
<^p/ <^u> — 1.15±0.15, consistent with values found in jet 
fragmentation The r]' yield is scaled to be consis- 

tent with jet fragmentation a^y /a,-i ~ 0.15 ± 0.1 5 [p3| . 
The ■0' is adjusted to be a^i jajj^ — 0.14 ± 0.03 [|68[ in 
agreement with PHENIX measurements p7| . 

For the ry, w, (/), and J/V", and also 77', and "0' (in P + p) 
the quoted uncertainties include those on the data as well 
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TABLE VI: 


Systematic uncertainties on tlie dilepton yield and mass range 


of applicability. 




Syst. Err. component 


P + P 


Au + An 


Mass Range 


Type 


pair reconstruction 


UA% 


13.4% 


0-4 GeV/c^ 


B 


conversion rejection 


6% 


6% 


0-0.6 GeV/c^ 


B 


pair cuts 


5% 


5% 


0.4-0.6 GeV/c^ 


B 


occupancy efficiency 


- 


3% 


0-4 GeV/c^ 


C 


BBC and trigger bias 


11.3% 


- 


0-4 GeV/c^ 


C 


ERT efficiency 


5% (20% for mr <1 GeV/c^) 




0-4 GeV/c^ 


B 


combinatoriai background 


3%-B/S 


0.25%-B/S 


0-4 GeV/c^ 


B 


correlated bacfcground 


mass-dependent (Fig. p^) 


mass-dependent (Fig. 


0-4 GeV/c^ 


B 


centrafity 




10% 


0-4 GeV/c^ 


B 


acceptance correction 


10% 


10% 


0-4 GeV/c^ 


B 


charm acceptance 


5% 


5% 


>0.5 GeV/c^ 


B 




(GeV/c) 



p.^ (GeV/c) 



FIG. 24: (Color online) Compilation of meson production cross sections in p -I- p (left) and Au + Au (right) collisions at 
y^SNN = 200 GeV. Shown for p + p are data for neutral and charged pions [psil , rj [ |56K kaoris [HsI , lo |^^, (j) S], and 
J/ip [^. Shown for Au -I- Au are data for neutral and charged pions [^l|, 77 [ |56| , kaons |3l|, oj, |p2|, and J/ip [ |63{ . The 
data are compared to the parameterization based on tot scaling used in our hadron decay generator. 



as those using different shapes of the px distributions 
to extrapolate to zero pT- Specifically we have fit the 
functional form given in Equation 
free and also with an exponential c 



30 with all parameters 
istribution in mr- For 
p nor in Au + Au, 
Au, the 



the p, which is not measured in p 
and T]' and ip' , which are not measured in Au 
uncertainty in the table represent the quadrature sum of 
the uncertainty of the cross section and the uncertainty 
relative to other mesons. 

All the mesons shown in Fig. can be described by 
the TOT-scaling parameterization of the pion spectrum. 
The fact that the ry's and the kaons follow the same ttit- 



scaling prediction over all pt appears to be due to the 
fact that the masses of the particles are almost the same. 
In Au -|- Au however rj and kaons do not follow the same 
TOT-scaling prediction. At high px, where we measure 
?7's, we see that they are suppressed as much as pions 
and the trend of Raa for these two mesons looks iden- 
tical. However we have observed a different trend (i.e. 
a smaller suppression) for strange particles and 77 has a 
strangeness content too [t^]. Therefore, since the 
77 cannot be measured at low-pr, we take as systematic 
uncertainty in the \ow-pT region (and consequently on 
the extrapolated dN/dy), the difference of the two spec- 
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TABLE VII: Fit parameters from the modified Hagedorn 
function (Equation wCj) for p + p and Au + Au pion spectra 
(tt" and tt*) and tlie corresponding rapidity density dN/dy. 



Parameter 



p + p 



Au + Au 



dN/dy 
A[mh GeV-^c^ 
a[(GeV/c)-'l 
b[(GeV/c)-^l 
Po[GeV/c] 



1.06±0.11 
377±60 
0.356±0.014 
0.068±0.019 
0.7±0.02 
8.25±0.04 



95.7±6.9 
504.5±10 
0.52±0.007 
0.16±0.010 
0.7±0.005 
8.27±0.02 



tra (~ 30%) in the low-pr range. We note that this is a 
conservative estimate of the systematic uncertainty. Sta- 
tistical models which reproduce well the particle spectra 
and ratios measured at RHIC ^ calculate a dN/dy 
for the 77 which is well within the assigned systematic 
uncertainty. At high-pT we assign a smaller systematic 
uncertainty of 7% (17%) for p + p (Au + Au) collisions 
arising from the asymptotic value of the ry/Tr*^ ratio of 
0.48 ± 0.03 (0.08) based on PHENIX measurement 

Once the meson yields and pt spectra are known 
the dilepton spectrum is given by decay kinematics and 
branching ratios, which are implemented in our de- 
cay generator EXODUS following earlier work published 
in [ p6[ . The branching ratios are taken from the compi- 
lation of particle properties in |^ . For the Dalitz decays 
TT*^, ?7, 77' — > e~'"e~7 and the decay lo —> e^e~Tr^ we use the 
Kroll-Wada expression |7^ with electromagnetic transi- 
tion form factors measured by the Lepton-G collabora- 
tion ItJ]. For the decays of the vector mesons p, cu, 
— > e+e~ we use the expression derived by Gounaris and 
Sakurai extending it to 2 GeV/c^, slightly beyond 

its validity range. For the J/V' and ip' e^e~ we use 
the same expression discussed in ||5^ modified to include 
radiative corrections. All vector mesons are assumed to 
be unpolarized. For the Dalitz decays, where the third 
body is a photon, the angular distribution is sampled ac- 
cording to 1 -I- cos^ dcs, where 0cs is the polar angle of 
the electrons in the Collins-Soper frame. 

The resulting systematic uncertainties on the mass 
spectrum depend on mass, and range from 10 to 30%. 
They result primarily from the uncertainty on the mea- 
sured pion yield and on the meson-to-pion ratios. The 
uncertainty from the measured electromagnetic transi- 



tion form factors, in particular for the w 



de- 



cay, is also included but contributes significantly only in 
the range around 0.5 to 0.6 GeV/c^. The uncertainty 
from polarization is negligible but also included. 



V. RESULTS 

This Section presents the results for the p + p and the 
Au -I- Au analyses. The p + p data provide a good base- 
line for understanding the results of the Au -I- Au analy- 



sis presented in this Section. This Section is organized as 
follows. Each subsection concentrates on a different re- 
gion of the (mee , Pt) phase-space. For each region we will 
present th e resu lts for the p+p and for the Au -I- Au data. 
In Section V A wc will show the inclusive mass spectrum 
for p + p and minimum bias Au -I- Au collisions and we 
will compare it with the cocktail. In Section |VB| we dis- 
cuss the results in the IMR by comparing the data with 
the charm exp ectations from (A'^coirscaled) pythia ^7\ . 
In Section V C we present the yields in the LMR and dis- 
cuss their centrality dependence. Section VD discusses 



the Pt dependence of the mass spectra. Part of the LMR, 
denoted with LMR I, where the pt of the e+e" pair is 
much larger than its mass, is the region of quasi-real vir- 
tual photons: we present the resu lts i n this re gion for 
p + p and Au -t- Au data in Section VE. Section VF dis- 



cusses what we can learn from the measurement of direct 
photons (in LMR I) about the yields in the pr-hiclusive 
LMR. Finally Section V G presents the pt spectra for dif- 
ferent mass bins for p + p and Min. Bias Au + Au data 
and compares them with the expectations from cocktail 
plus charm plus direct photons. 



A. pT-inclusive Mass Spectra 



Figure 25 compares the yield of e+e" pairs in the 
PHENIX acceptance in p + p data to the expected yield 
from the contributions of the cocktail and shows the var- 
ious sources of the cocktail (hadron decays, charm, bot- 
tom and Drell-Yan pairs) . The p + p data are very well 
described by the expectation from the hadronic cocktail 
and heavy flavor decays for the entire mass range within 
the uncertainty of the data and the cocktail. 

Figure ^ compares the e+e" yield in the PHENIX 
acceptance in minimum bias Au -I- Au collisions to the 
expected yield from the contributions of various sources. 
The cocktail sources are the same as in the p + p data, 
but tuned separately to the Au + Au measurements. 
The data below the pion mass, where 7r° Dalitz decays 
dominate, are well described by the cocktail. The vec- 
tor meson yields in the cocktail are partly based on the 



e+e pair data (see Sec. IV) and consequently agree well 



with the data. In particular the J/tp yield in the cock- 
tail is exclusi vely based on the e+e^ measurement [6^ 
(shown in Fig, 24 and therefore the data in Fig. 26 are ex- 



pected to agree with the cocktail at the J/tp peak. J/^p 
suppression would be observed if the cocktail would in- 
stead use the A^coirscaled p + p measurement js^ . In 
the IMR, besides the charm, bottom and Drell Yan cal- 
culated with PYTHIA which are the same as in the 
p + p data, scaled by A^coU, there is another curve drawn 
for the charm. This will be described in Section |VB| . 

Figure ^ shows the mass spectra for p + p (bottom) , 
for minimum bias Au -I- Au (top) and for five differ- 
ent centrality classes in Au + Au. The data are com- 
pared to the sum of hadronic cocktail and charmed me- 
son decays. The charm cross section, measured in p + p 



29 



TABLE VIII: Hadron rapidity densities used in our hadron decay generator. For the uj and 
used together with data from the quoted references. 



dN 1 



relative uncertainty. 



mcson/vr" 



data from this analysis were 



data used 



P 

v' 



1.065 ±0.11 
(1.1 ±0.3) X 10-^ 
(8.9 ±2.5) X 10"^ 
(7.8 ±1.8) X 10^2 
(9.0 ±2.0) X 10~^ 
(1.3 ±0.5) X 10^^ 
(1.77 ±0.27) X 10"^ 
(2.5 ±0.7) X lO"** 



P + P 
10% 
30% 
28% 
23% 
24% 
40% 
15% 
27% 



1.0 

1.032 X 10"^ 
8.34 X 10"^ 
7.32 X 10"^ 

8.4 X 10"^ 
1.27 X 10"^ 
1.66 X lO"'' 

2.3 X 10"® 



PHENIX |54j5|] 
PHENIX|||] 
jet fragmentation ||53 
PHENIX |57 
PHENIX 158 
PHENIX 
PHENIX 15C 

PHENIX 




IT 

V 

p 

UJ 

V' 

J/i> 



(9.572 ±0.95) X 10 
(1.077 ±0.32) X 10 

8.60 ±2.8 

9.88 ± 3.0 

2.05 ±0.6 

2.05 ±0.2 
(1.79 ±0.26) X 10"^ 
(2.6 ±0.7) X 10"* 



Au ± Au 
10% 
30% 
33% 
30% 
30% 
100% 
15% 
27% 



1.0 

1.12 X 10"* 
8.98 X 10"^ 
1.03 X 10"* 

2.14 X 10"^ 

2.15 X 10"^ 
1.82 X lO"'^ 
2.70 X lO"'' 




PHENIX 
PHENIX |5 
jet fragmentation 
PHENIX Q 

PHENIX jell 

PHENIX and 
PHENIX Hi 
PHENIX and 



Uc-c = 567±57'^'^*±224"y''' ||l has been scaled by iVcou 
(given in Table |). For each centrality class, the data and 
the cocktail are absolutely normalized. Each data set is 
compared with two corresponding cocktail lines, shown 
in solid and dotted curves. The difference between the 
cocktails is due to uncertainty in the cc contribution (see 
discussion below). 

Unlike the p + p mass spectrum, the Au + Au mass 
spectra show enhancement above the cocktail, in partic- 
ular for the LMR (0.15 - 0.75 GeV/c^). There is lit- 
tle enhancement for peripheral (60-92%) data, but very 
strong enhancement for two most central classes (0-10% 
and 10-20%). The enhancement increases rapidly with 
increasing centrality. 

In order to quantitatively describe this enhancement, 
more information is needed about other components that 
can potentially contribute to the LMR, namely the open 
heavy flavor and internal conversion of real direct pho- 
tons. We discuss them in the next sections. 



B. Open Heavy Flavor Contribution 

The dilepton yield in the IMR is dominated by semi- 
leptonic decays of charm hadrons correlated through fla- 
vor conservation. Small contributions also arise from bot- 
tom hadrons and Drell-Yan. For p-\-p data we determine 
the heavy flavor contribution by subtracting the hadronic 
cocktail from the dilepton data. We integrate the sub- 



tracted yield in the IMR, extrapolate to zero e+e" pair 
mass to get the entire cross section, correct for geometric 
acceptance, and convert to a production cross section us- 
ing known branching ratios of semi-leptonic decays [ p3| . 
Details of the analysis of the charm cross section are re- 
ported in 

We find a rapidity density of cc pairs at mid rapidity: 

^|„=o = 118.1 ± 8.4"*"* ± 30.7 syst ± 39.5'"°^"Vb. 

dy 

This corresponds to a total charm cross section of CTcc — 
544±39"*''*± 142">"'' ±200'"°'^'=Vb, consistent with previ- 
ous measurement of single electrons by PHENIX (ctcc = 
567±57"*"'±224''y"* ^b) @ and with a fixed-order-plus- 
next-to-leading-log (FONLL) pQCD calculation (dee — 
256l«°Mb) |7§. 

In Au ± Au the dynamic correlation of c and c, which is 
essential to determine the mass spectral shape, could be 
modified compared to p±p collisions. The observed sup- 
pression and the elliptic flow of non-photonic electrons in- 
dicates that charm quarks interact with the medium 
which should change the correlations between the pro- 
duced cc pairs. We also note that the pT distribution 
for electrons generated by pythia is softer than the 
spectrum measured lap + p data but coincides with that 
observed in minimum bias Au ± Au collisions. Thus we 
compare our Au ± Au data to two extreme scenarios that 
bracket the charm contribution: 

(i) The correlation is unchanged by the medium and 
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FIG. 25: (Color online) Inclusive mass spectrum of e^e~ 
pairs in the PHENIX acceptance in p + p collisions compared 
to the expectations from the decays of light hadrons and cor- 
related decays of charm, bottom, and Drell-Yan. The contri- 
bution from hadron decays is independently normalized based 
on meson measurements in PHENIX. The bottom panel shows 
the ratio of data to the cocktail of known sources. The sys- 
tematic uncertainties of the data are shown as boxes, while 
the uncertainty on the cocktail is shown as band around 1. 



equals what is known from p-fp collisions. In this case we 
can use the same pythia calculation scaled to match the 
cross section measured 'm p + p and scale it by the mean 
number of binary N + N collisions (as given in Table |). 

(ii) The cc dynamical correlation is washed out by 
medium interactions, i.e. the direction of c and c quarks 
are uncorrelated. We sample from the heavy flavor 
single-electron pT spectra, choose the angle randomly 
and keep the overall cross section fixed to the experi- 
mental data 1^. Because the average opening angle of 
uncorrelated pairs is smaller than the one resulting from 
the back-to-back correlation predicted by PYTHiA, the 
mass spectral shape of uncorrelated pairs is much softer 
than the one calculated by pythia. 

The charm contribution determined by case (i) is 
shown as the upper dashed curve in Fig. ^ or the upper 
solid curves in Fig. ^ The charm contribution deter- 
mined by case (ii) is shown as the dotted curve in Fig. 
or the dashed curves in Fig. |2^. In both cases the total 
yield of charm is normalized to the value measured by 
PHENIX. 

In the Min. Bias Au + Au data set, the IMR seems to 
be well described by the continuum calculation based on 
case (i). This is somewhat surprising, since single elec- 
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FIG. 26: (Color online) Inclusive mass spectrum of e+e~ 
pairs in the PHENIX acceptance in minimum-bias Au + Au 
compared to expectations from the decays of light hadrons 
and correlated decays of charm, bottom and Drell-Yan. The 
charm contribution expected if the dynamic correlation of c 
and c is removed is shown separately. Statistical (bars) and 
systematic (boxes) uncertainties are shown separately. The 
contribution from hadron decays is independently normalized 
based on meson measurements in PHENIX. The bottom panel 
shows the ratio of data to the cocktail of known sources. The 
systematic uncertainties of the data are shown as boxes, while 
the uncertainty on the cocktail is shown as band around 1. 

tron distributions from charm show substantial medium 
modifications Thus, it is hard to understand how the 
dynamical correlation at production of the cc remains un- 
affected by the medium. Case (ii) leads to a much softer 
mass spectrum, as shown by the dotted curve in Fig. p6| ). 
This would leave significant room for other contributions, 
e.g. thermal radiation. 

We have integrated the yield in the mass region 1.2 to 
2.8 GeV/c^ and normalized to the number of binary col- 
lisions iVcoii (Fig. ^8|). The systematic uncertainty due 
to A^coU (as indicated in Table |) has been included in 
the overall systematic uncertainty. Within uncertain- 
ties iVcoii scaling is observed for the production of non- 
photonic electrons, i.e., for those electrons arising from 
decays of heavy-flavor hadrons |6| . The normalized yield 
shows no significant centrality dependence and is consis- 
tent within systematic uncertainties with the expecta- 
tion based on A^coii-scaled pythia, with the cross section 
measurement of (case (i)). However the scaling with 
A'coii may be a mere coincidence resulting from two bal- 
ancing effects: the energy loss of charm, which increases 
with iVpart, would lead to a softer mass distribution and 
therefore less yield in the IMR (case (ii)), while a ther- 
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p+p & Au+Au at\/ifjM = 200 GeV 

• Au+Au min. bias x 10^ 

■ Au+Au 0-10% X 10^ 

• Au+Au 10-20% X 10 

■ Au+Au 20-40% 

• Au+Au 40-60% X 10"' 
Au+Au 60-92% X 10"^ 
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FIG. 27: (Color online) Invariant mass spectrum of e+e^ 
pairs inclusive in pr compared to expectations from the model 
of hadron decays for p + p and for different Au + Au central- 
ity classes. The charmed meson decay contribution based on 
PYTHIA j??! is included in the sum of sources (solid black 
line). The uncorrelated charm is shown with the dotted 
line. Statistical (bars) and systematic (boxes) uncertainties 
are shown separately. The systematic uncertainty on the ex- 
pected hadronic sources is not shown: it ranges from ~10% 
in the tt" region to ~30% in the region of the vector mesons. 
The uncertainty on the charm cross section, which dominates 
the IMR, is ~30% both in p -f p and in Au + Au collisions. 



mal contribution could increase faster than linearly with 
-^part resulting in more yield in the IMR. Such a coinci- 
dence may have been observed at the SPS Q, where 
a prompt component has now been suggested by NA60 



C. Low-Mass Excess in Au + Au Data 



Figure 27 shows that the low-mass enhancement is con- 
centrated in the first two centrality classes, i.e. 0-10% 
and 10-20%. For more peripheral collisions the enhance- 
ment diminishes. Only some small excess is visible for 
20-40% and 60-92% while no deviation is observed in 
40-60% with respect to the cocktail beyond systematic 
uncertainties. 

To quantify the centrality dependence of the enhance- 
ment, we have integrated the yield in two mass win- 
dows: below 0.1 GeV/c^ and 0.15 to 0.75 GeV/c^ Since 
the cocktail yield in these regions arises mostly from 
hadrons (more than 90% from tt" below 0.1 GeV/c^ 
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FIG. 28: (Color online) Dielectron yield per binary collision 
in the mass range 1.2 to 2.8 GeV/c^ as a function of A'part. 
Statistical and systematic uncertainties are shown separately. 
Also shown are two bands corresponding to different estimates 
of the contribution from charmed meson decays. The width of 
the bands reflects the uncertainty of the charm cross section 
only. 
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FIG. 29: (Color online) Dielectron yield per participating 
nucleon pair (Afpart/2) as function of A^part for two different 
mass ranges (a: 0.15 < ruee < 0.75 GeV/c^, b: < mee < 
0.1 GeV/c^) compared to the expected yield from the hadron 
decay model. The two lines give the systematic uncertainty 
of the yield from cocktail and charmed hadron decays. For 
the data statistical and systematic uncertainties are shown 
separately. 



and more than 99% from the sum of light hadrons in 
0.15 < TUee < 0.75 GeV/c^) we compare the measured 
yield to the rate of pion production. Pions were found 
to scale approximately with iVpart |7q] , therefore we com- 
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pare the measured yield in data to A^part- 

The top panel of Fig. ^ shows the centrality depen- 
dence of the yield in the mass region 0.15-0.75 GeV/c^ 
divided by the number of participating nucleon pairs 
(iVpart/2). The systematic uncertainty due to iVpart (as 
indicated in Table |) has been included in the overall 
systematic uncertainty. For comparison the yield below 
0.1 GeV/c^, which is dominated by low-p^ pion decays, 
is shown in the lower panel. 

For both mass intervals the yield is compared to the 
yield calculated from the hadron cocktail. Two solid 
curves on each panel show the upper and lower limit of 
the expected yield from the cocktail. The cocktail uncer- 
tainty includes the uncertainty in the charm contribution 
discussed in the previous section. In the lower mass range 
the yield agrees with expectations, and is proportional to 
the pion yield (bottom panel of Fig. |2^). In contrast, in 
the range from 0.15 to 0.75 GeV/c^ the observed yield 
rises significantly above expectations. 

The enhancement factor, defined as the ratio between 
the measured yield and the expected yield for 0.15 < 
ruee < 0.75 GeV/c^, is 4.7 ± 0.4''*'^* ± 1.5"^"* ± 0.9'°°'^°', 
for Min. Bias data. The first error is the statistical error, 
the second the systematic uncertainty of the data, and 
the last error is an estimate of the uncertainty in the 
cocktail, i.e. the expected yield from hadronic sources. 
For the various centrality bins the enhancement factor is 
reported in Table 
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TABLE IX: The enhancement factor, defined as the ra- 
tio between the measured yield and the expected yield for 
0.15 < rriee < 0.75 GeV/c^, for different centrality bins. The 
meaning of the errors is defined in the text. 

Centrality Enhancement (zbstat ±syst zbmodel) 



00-10 % 
10-20 % 
20-40 % 
40-60 % 
60-92 % 
Min. Bias 



7.6 ± 0.5 ± 1.5 ± 1.5 
3.2 ± 0.4 ± 0.1 ± 0.6 

1.4 ± 1.3 ± 0.02 ± 0.3 
0.8 ± 0.3 ± 0.03 ± 0.2 

1.5 ± 0.3 ± 0.001 ± 0.3 

4.7 ± 0.4 ± 1.5 ± 0.9 



The increase is qualitatively consistent with the con- 
jecture that an in- medium enhancement of the dielectron 
continuum yield arises from scattering processes like tttt 
or qq annihilation. In this case the enhancement would 
scale proportional to N^^^.^^, different from the hadronic 
cocktail which scales proportional to A'part- 



D. pt Dependence of the Mass Spectra 

Figure |30| compares e+e^ invariant mass spectra mea- 
sured in p+p and in Min. Bias Au -I- Au collisions to the 
corresponding expectations from the cocktail of hadron 
decays and open charm, in different ranges oipx- Data 
and cocktail are absolutely normalized. 



The "knee" in the mass spectra beginning at rUee — 
0.1 GeV/c^ corresponds to the cut-off of the 7r° Dalitz 
decays. For rriee > ni^, the 77 Dalitz decay (77 — > e~^e~j) 
is the dominant hadronic source of e"'"e~ pairs, followed 
by the u Dalitz decay (a; — > e"'"e~7r''), with small contri- 
butions from other sources such as rj' and (f>. The detector 
acceptance and resolution effects broaden the low-mass 
Dalitz peak and smear the "knee" of the 7r° Dalitz contri- 
bution. These detector effects are included in the cocktail 
calculation. 

The p + p data are consistent with expectations from 
the cocktail over the full mass range in the low px bin. In 
the highest px bins however the data are enhanced with 
respect to the cocktail. The deviations are however small 
in contrast to the Au -t- Au data which show a large en- 
hancement in the LMR above , which is concentrated 
at low-pT- For px > 1.0 GeV/c the enhancement be- 
comes smaller than at lower px but it is still larger than 
the one observed in the p + p data. 

In the introduction, we classified the low mass region 
into LMR I (low mass, high px) and LMR II (low mass, 
low Px) (see Fig. |^). The behavior of the low mass excess 
in Au -I- Au shown in Fig. ^ is different for LMR I and 
LMR II. In LMR I, the excess has a similar shape to 
the cocktail and the level of the excess with respect to 
the cocktail is approximately constant. In this region, 
a contribution from internal conversion of virtual direct 
photons is expected. In LMR II, the excess increases with 
increasing mass and decreasing px ■ 

In the following we first analyze the data in the low 
mass high px region (LMR I), which allows the measure- 
ment of direct photons; then we study how much they 
contribute to the inclusive enhancement, dominated by 
the yield in the LMR II. 



E. Measurement of Direct Photons 

In LMR I (Fig. |^), where the px of the e+e" pair is 
much greater than its mass (mee ^ Pt), the yield of 
the virtual photons is approximately the same as that of 
real photons. Therefore, in this quasi-real virtual pho- 
ton region, the production of real direct photons can 
be deduced from measurements of e~^e~ pairs. Theo- 
retical details are given in Appendix ^ In this Section, 
we determine the real direct photon cross section for 
1 < Pt < 5 GeV/c from the data shown in Fig. |3^. 
We use the mass range 0.1 < iriee < 0.3 GeV/c^. 

The relation between real photon production and the 
asso ciat ed e+e~ production can be written as (Equa- 
tion IbI) 



dmei.dpx 37r m 



2a 1 ^dN^ 



L(mee) = W 1 5- 



Here a is the fine structure constant, 



(31) 

(32) 
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FIG. 30: (Color online) e^e~ pair invariant mass distributions in p + p (left) and minimum bias Au + An collisions (right). 
The pt ranges are shown in the legend. The solid curves represent the cocktail of hadronic sources (see Section IV ) and include 
contribution from charm calculated by pythia using the cross section from scaled by A^coii. 
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FIG. 31: (Color online) Electron pair mass distribution for 
Au + Au (Min. Bias) events for 1.0 < pr < 1.5 GeV/c. 
The two-component fit is explained in the text. The fit range 
is 0.12 < < 0.3 GeV/c^. The dashed (black) curve at 
greater m^e shows f{mee) outside of the fit range. 



of e+e^ pair, me is the electron mass, and S{mee,PT) 
is a process-dependent factor that accounts for differ- 
ences between real and virtual photon production, such 
as form factors, phase space, and the spectral function. 
Equation ^ holds for any process emitting real photons, 
in particular direct or thermal emission. For high px 
(Pt ^ iTiee) the process dependence becomes negligi- 
ble and the factor S{mee,PT) becomes 1 as iriee 



or rriee/pT 0. For rriee 3> the factor Liniee) also 
becomes very close to unity. Thus the relation simplifies 
to 



2a 1 dN^ 



drrieedpT 37r TOee dpT 



(33) 



Here the mass distribution of electron pairs for a given pT 
bin takes on a very characteristic 1 / niee shape. If there 
is real direct photon production in a given px bin, there 
should be a corresponding electron pair contribution that 
behaves like l/ruee in the same pt bin. Therefore, the 
real direct photon production can be determined from 
the yield of the excess electron pairs. 

For Dalitz decays, the l/niee behavior is truncated by 
the kinematic limit and S{mee) becomes zero for rriee > 
mil, where rrih is the mass of the hadron. The functional 
form of SirUee) for Dalitz decays is given in Appendix^. 
In contrast, the factor S{mee,q) for the direct photon 
process is unity for pT ^ Wee . We exploit this difference 
to separate the direct photon signal from the hadronic 
background. Since 80% of the hadronic photons are from 
tt'^ Dalitz decays, the signal to background (S/B) ratio 
for the direct photon signal improves by a factor of five 
for TUee > m^o « 0.135 GeV/c^, thereby allowing a real 
direct photon signal that is 10% of the yield of hadronic 
decay photons to be observed as a 50% excess of e+e~ 
pairs for this mass range. 

Figure ^ shows a visible excess above the tt" cutoff for 
all Pt bins of the Au + Au data. For pT > 1 GeV/c, the 
excess is almost a constant factor above the cocktail. As 
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we examine later the mass distribution for pT > 1 GeV/c 
is consistent with the 1 /niee shape expected for the elec- 
tron pairs from internal conversion of virtual direct pho- 
tons. 

In the following we assume the excess for pT > I GcV/c 
and niee < 0.3 GeV/c^ is entirely due to internal conver- 
sion of virtual direct photons and deduce the real di- 
rect photon yield from the e+e" pair yield using Equa- 
tion We demonstrate the validity of this assump- 
tion later. Although the data are consistent with l/rngg 
over a wider mass range (mee ~ 0.7 GeV/c^), as can 
be seen in Fig. we limit our analysis for 0.1 < rriee < 
0.3 GeV/c^. We do so in order (i) to ensure the condition 
"lee <C PT for the lowest pr bin (1.0 < pr < 1-5 GeV/c), 
(ii) to keep the correction factor S{M, q) close to unity, 
and (iii) to minimize uncertainty due to cc. In this kine- 
matic range, the contribution of cc e+e^, estimated 
by PYTHIA, is less than 5% of the excess. 

In order to quantify the excess, we fit a two-component 
function, 

/(mee;r) = (1 - r)/c(mee) -I- r/dir(mee) (34) 

to the mass distribution. Here /c("iee) is the shape of the 
cocktail mass distribution (shown in Fig. /dir('7iee) is 
the expected shape of the virtual direct photon internal 
conversion mass distribution, and r is the only fit parame- 
ter. In the low mass region used for the fit, the functional 
form of fcijnee) is the sum o f D alitz decay rnass distribu- 
tions of hadrons (Equation BS - Equation |B10| ) filtered 
through the PHENIX acceptance and smeared by the 
detector effects. It is calculated by a Monte Carlo sim- 
ulation which takes into account detector effects such as 
finite mass resolution. The functional form of /dir(™ee) 
corresponds to Equation ^ with S'(mee) = 1. It is also 
filtered through the PHENIX acceptance and smeared by 
detector effects. 

Both fc{mee) and /dir("iee) are separately normalized 
to the data for mee < 30 MeV/c^. In this mass region 
Sirriee) of 7r° Dalitz decays is very close to unity. Thus 
the functional shapes of /c and /dir are essentially iden- 
tical and equal to L{mee)/'m.ee smeared by the detector 
effects. This means that the fit function f{mee',r) in 
this mass range is independent of the fit parameter r as 

(l-r)L(mee)/mee-|-rL(mee)/TOee = L{mee) / rUee- ThuS 

this normalization ensures that the yield of fit function 



is always normalized to that of the data for 



TOee < 30 MeV/c^. The parameter r can be interpreted 
as the direct photon fraction of the inclusive photon yield. 

This fitting method has the advantage of canceling 
most of the systematic uncertainties of the cocktail nor- 
malization relative to the data. The PHENIX accep- 
tance for electron pairs with uriee < 0.3 GeV/c^ and 
Pt > I GeV/c is almost constant, and its shape can be 
calculated accurately as a function of mass. Many sys- 
tematic effects, such as electron identification efficiency, 
detector dead area, etc. can influence the absolute value 
of the acceptance but not its shape. 

For each px bin, /(mee) is fit to the data for several 



mass ranges with r the only fit parameter. Figure |3l| 
shows /dii('7iee) and fcipiee) together with the fit result 
for Au -I- An Min. Bias data for 1.0 < < 1-5 GeV/c 
and the cocktail components. The dashed curve shows 
/("lee) extended outside of the fit range. Although the 
mass region is not used in the fit, the fit function de- 
scribes the data for TOee > 0.3 GeV/c^. 

The fit shown in Fig. |l| has xVNDF= 12.2/6. The 
somewhat large values is due to the large contribu- 
tion from the lowest mass bins, where statistical errors 
are small and systematic errors due to the detector res- 
olution are significant. The value is calculated from 
the statistical errors only. The results for the fit range 
0.12 < TOee < 0.3 GeV/c^ are summarized in Table 
For pt > 1.5 GeV/c^ the fit gives good x^/NDF, demon- 
strating that the shape of the excess is consistent with 
1 /TOee as expected for internal conversion. 

To evaluate the systematic uncertainty due to the mass 
range used for the fit, the fit was repeated for three 
mass ranges: 0.08 < TOee < 0.3 GeV/c^, 0.1 < TOee < 
0.3 GeV/c^, and 0.12 < TOee < 0.3 GeV/c^. The value of 
r is taken as the average of the results for these three fit 
ranges. 

The sources of systematic uncertainty on the fit in- 
clude (1) the fit range, (2) the mass spectrum of the 
data, and (3) the cocktail. The sources of the system- 
atic uncertainty on the mass spectrum relative to the 
cocktail include (1) uncertainties on the correlated back- 
ground due to jet pairs (~ 2% for p + p), (2) uncertain- 
ties in the acceptance and efficiency in the mass range 
of the signal (0.1 < mee < 0.3 GeV/c^) relative to the 
mass range (mee < 30 MeV/c^) used for the normal- 
ization (~ 1% for acceptance and ~ 1% for efficiency), 
and (3) uncertainty in the mixed-event normalization 
{0.25%/ (S/B)), where S/B is the signal-to-background 
ratio in 0.1 < TOee < 0.3 GeV/c^. The uncertainty (2) 
is small since uncertainties in the absolute normalization 
cancel when the cocktail is normalized to the data in the 
low- mass peak (TOee < 30 MeV/c^). The largest source 
of uncertainty is the particle composition in the hadronic 
cocktail, namely the rj/n'^ ratio. This corresponds to a 
~ 7% (~ 17%) uncertainty in the p+p (Au -I- Au) cocktail 
for 0.1 < TOee < 0.3 GeV/c^. All systematic uncertainties 
are added in quadrature to obtain the total systematic 
uncertainty. 



TABLE X: Summary of the fits to Equation |34| in the range 
0.12 < mee < 0.3 GeV/c^. 



PT(GeV/c) 



xVndf 



1.0-1.5 
1.5-2.0 
2.0-2.5 
2.5-3.0 
3.0-4.0 
4.0-5.0 



0.189 ±0.021 
0.165 ±0.022 
0.146 ±0.029 
0.165 ±0.040 
0.224 ± 0.048 
0.206 ± 0.093 



12.2/6 
4.6/6 
6.6/6 
3.3/6 
3.7/6 
4.2/3 
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Since the r?/7r° ratio is the largest source of uncertainty, 
we also studied fits with a three component function, 

hirriee) = (1 - r - ry^)fc{mee) + rfdh {mee) + rnfnimee), 
with a constraint on r,, such that rj/ir'^ — 0.48±0.03(0.08) 
for p + p (Au + Au) ll^. These alternative fits give 
consistent results for r within statistical uncertainties. 

So far we have assumed that the excess in pT > 
1 GeV/c and 0.1 < Wee < 0.3 GeV/c^ is entirely due to 
internal conversion of virtual direct photons. This means 
that we assume S{mg^g, q) for excess virtual photons to be 
unity. In the following, we examine the validity of the as- 
sumption from the data. 

As shown in Appendix the shape of the virtual pho- 
ton spectrum as a function of mass can be obtained from 
the electron pair yield as 



dN* 37r drill 



(35) 



Since the shape of /dir(TOee) is i/rUee smeared by the 
detector effects, a fit of i? = (data — cocktail) //dir(TOee) 
to a constant can be used to test that the excess has the 
shape expected for internal conversion of direct photons. 
Note that the detector effects in the numerator and de- 
nominator cancel in the ratio. 

Furthermore, since /dir(TOee) is normalized to the data 
for rriee < 30 MeV/c^, R can be interpreted as the ratio 
of the virtual photon yield to the inclusive real photon 
yield: 



R{m,pT) 



(to,Pt) dN"'^\pT) 



dpi 



-I 



dpi 



(36) 



5(m,pT)diVf ^^*(PT)/diVi"^'(PT).(37) 



Figure |3^ shows i? as a function of rriee for 0.1 < 
< 0.7 GeV/c^. The ratio cannot be measured for 
rriee < 0.1 GeV/c^ because in this mass region the signal 
is masked by large background from tt" Dalitz decays. 
The distributions are consistent with a constant for the 
five pt bins. For the highest pT bin (j)t > 4 GeV/c) 
the shape of the virtual photon mass spectrum is not 
well constrained due to limited statistics. However, it 
is reasonable to expect that the same constant behavior 
continues for higher pT- The ±lcr band of a constant 
value fit are shown in each panel. Table XI summarizes 
the results of the fits. For all pr bins, the constant value 
fit gives a good x^/NDF value. This demonstrate that 
the data are consistent with a constant 5 (mee) for these 
Pt bins. 

As discussed in Appendix the ratio is expected 
to be a smooth function of rriee- Hadronic and par- 
tonic direct photon contributions to S {rriee, q) are ex- 
pected to be nearly constant in this range, qq annihi- 
lation can make a contribution proportional to Tn^g, but 
should be much smaller than these two components for 
0.1 < mee < 0.3 GeV/c^ @ |ll |§. There is no sign of 
a component that scales with m^^ for rriee < 0.3 GeV/c^, 
suggesting that the qq contribution is indeed small. Fig- 
ure a2 illustrates that S{M, q) is indeed constant, and 
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FIG. 32: (Color online) Ratio R = 

(data — cocktail) //dir(»Tiee) of electron pairs for different 
Pt bins in Min. Bias Au -I- Au collisions. The pt range of 
each panel is indicated in the figure. The yellow band in each 
panel shows ±1(t range of a constant fit value to the data 
points. 



supports the use of r obtained from our fit as the direct 
photon fraction in the inclusive photon spectrum. 

The absence of any increase in R{mee) for rriee > 
0.5 GeV/c^ is somewhat surprising. A constant R as 
function of ruee implies that the S{mee) is also a con- 
stant. If the excess electron pairs are thermal pairs from 
the medium, we expect an increasing contribution from 
qq annihilation or from the tail of the (possibly modified) 
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FIG. 33: (Color online) The fraction of the direct photon 
component as a function of pr- The error bars and the error 
band represent the statistical and systematic uncertainties, 
respectively. The curves are from a NLO pQCD calculation 
(see text). 



p resonance, which leads to increase of Siniee) in higher 
mass. The data show no indication of such increase. For 
thermal radiation, S{mee) is the space time average of 
a product of electromagnetic sp ectra l function and the 
Boltzmann factor, see Equation B24 . We note that the 
Boltzmann factor — l/(e^/-^ — 1) can cause significant 
suppression of S{niee, q) for niee > 0.5 GeV/c^. Thus the 
nearly constant behavior of R{mei,) for m^ee > 0.3 GeV/c 
should not be interpreted as absence of any contribution 
other than internal conversion of virtual direct photons 
at large rrie 
Figures p 

of the direct photon component determined by the two- 
component fit (see Equation |3^) in p + p and Au + Au 
collisions, respectively. The curves represent the expec- 
tations from a next-to- leading-order perturbative QCD 
(NLO pQCD) calculation g^. For p + p, the curves 
show the ratio da^ / da'^''\pT) , and da'^^^ipT) is 
the p + p inclusive photon cross section (obtained from 
the data as described later). For Au + Au, the curves 



I (a) and (b) show the fraction r = . '^'r^* ^ 

I ^ ^ ' inclusive 7 



represent rAArfcr;^ (^'T)/dA^" (pt), where Taa is the 
Glauber nuclear overlap function and dN^^'^^{pT) is the 
Au -f- Au inclusive photon yield. The three curves corre- 
sponding (from top to bottom) to the theoretical scales 



TABLE XI: Summary of a constant fit to the ratio data shown 
in Fig. lij The fit range is 0.11 < rUee < 0.7 GeV/c^ 



PT(GeV/c) 



{R) 



1.0 

1.5 ■ 
2.0 ■ 
2.5 ■ 
3.0 



1.5 
2.0 
2.5 
3.0 
4.0 



0.173 ± 0.015 
0.149 ± 0.016 
0.125 ± 0.020 
0.115 ± 0.028 
0.210 ± 0.033 



16.9/12 
14.8/12 
14.7/12 
9.1/12 
4.4/12 



AuAu Min. Bias x10 
AuAuO-20%x10^ 
AuAu 20-40% x10 
p+p 
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FIG. 34: (Color online) Invariant cross section {p + p) and 
invariant yield (Au -I- Au) of direct photons as a function of 
Pt- The filled points are from this analysis and open points 
are from |7^, ^ . The three curves on the p+p data repre- 
sent NLO pQCD calculations, and the dashed curves show a 
modified power-law fit to the p + p data, scaled by Taa. The 
dashed (black) curves are exponential plus the Taa scaled 
p + p fit. The dotted (red) curve near the 0-20% centrality 
data is a theory calculation [M . 



set to /i = 0.5 Pt, Pt, and 2 pT, respectively, show the 
scale dependence of the calculations. While the frac- 
tion r is consistent with the NLO pQCD calculation in 
p + p, it is larger than the calculation in Au + Au for 
PT < 4.5 GeV/c. 

The direct photon fraction r in Fig. |3^ is converted to 
the direct photon yield using dN^"{pT) = rx d7V™'^'(pT)- 
Here we determine the inclusive photon yield for each pT 
bin from the yield of low mass e'^e^ pairs in the range 
niee < 30 MeV/c^ using the following method. The 
differential yield of electron pairs is related to that of 
photons by Equation The process dependent factor 
S(rnee,q) is unity within a few percent for any source of 
photon for mee < 30 MeV/c^. Thus the measured yield 
of electron pairs {N^^*^^) in niee < 30 MeV/c^ for a given 
Pt bin is proportional to that of inclusive photons in the 
same pt bin. 
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(38) 
(39) 
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Here e'^'^'^ represents the acceptance of PHENIX. The 
same relation holds for the cocktail calculation of pho- 
ton and electron pairs. 
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dpT 



(40) 



Here jV™'^'^*'''' is the yield of electron pairs for niee < 
30 MeV/c^ in the hadronic cocktail calculation, and 
^^cocktaiy^p^ is the yield of photons in the same cal- 
culation. Thus we have 
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The systematic uncertainty in the inclusive photon spec- 
tra equals the systematic uncertainty in N^^^'^, which 
is summarized in Table VI. The total systematic uncer- 
tainty in Nee is approximately 20%. 

In Fig. 34 the direct photon spectra thus obtained are 
compared with the direct photon data from Q and 
NLO pQCD calculations. The systematic uncertainty of 
the inclusive photon yield is added in quadrature with the 
systematic uncertainties of the data. The p + p data are 
shown as an invariant cross section using da — a^^^^dN 
using 4'^°^ = 42 mb. 

The direct photon data of this analysis are obtained 
from the yield of e^e~ pairs using Equation |3l| under the 
assumption S{mee,q) = 1 for 0.1 < niee < 0.3 GeV/c^. 
Although we have shown that our data are consistent 
with this assumption, what we actually measure is the 
yield of e'^e" pairs in this mass range. For completeness, 
we give the relation between the direct photon yield de- 
duced by the analysis and the electron pair yield that is 
actually measured. The relation between the real photon 
yield and th e e~^ e~ pair yield in this mass range is given 
by Equation B7 . Thus the yield of the excess e~^e~ pairs 
for 0.1 < TUee < 0.3 GeV/c^ can be obtained by multiply- 
ing the photon yield by a factor of |^ log ^ ~ 1.7x 10^^. 

The pQCD calculation is consistent with the p+p data 
within the theoretical uncertainties for pT > 2 GeV/c. A 
similarly good agreement is observed for n'^ |8^. The 
p + p data can be well described by a modified power- 
law function {App{l + p'^/b)~") as shown by the dashed 
curve in Fig. The Au -I- Au data are above the p + p 
fit curve scaled by Taa for p-r < 2.5 GeV/c, indicating 
that the direct photon yield in the low-p^" range increases 
faster than the binary-scaled p + p cross section. 

We fit an exponential plus the TAA-scaled p + p iit 
function (Ae^P^/^ + Taa x App{l + p'^/b)-'^) to the 
Au -|~ Au data. The only free parameters in the fit are 
A and the inverse slope T of the exponential term. The 
systematic uncertainties in T are estimated by changing 
the p + p &t component and the Au -I- Au data points 
within the systematic uncertainties. The results of the 
fits are summarized in Table XH, where A is converted 
to dN/dy for > 1 GeV/c. For central collisions, 
T = 221 ± l%stat.) ± l^{syst.)MeV. If an unmodified 
power-law function (cx p^") is used to fit the p + p spec- 
trum, we find n = 5.40 ± 0.15, and T = 240 ± 21 MeV. 



F. pt Dependence of Low Mass Excess 



min. bias Au+Au 

• 0.4<p^<0.6GeV/cx 10 
■ 0.6<p.^<0.8 GeV/c 

• 0.8<p^<1.0 GeV/cx 10^ 
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FIG. 35: (Color online) The e^e~ pair invariant mass distri- 
butions in minimum bias Au + Au collisions for the low-pT 
range. The solid curves represent the cocktail of hadronic 
sources (see Section IV ) and include contribution from charm 
calculated by pythia using the cross section from scaled 
by A^coii • 

The shape of the enhancement in Au -I- Au data in 
LMR II (low mass low px) is quite different from that in 
LMR I (low mass high px), where it behaves like l/ruee 
and is consistent with internal conversion of direct pho- 
tons. In LMR II, the enhancement is larger, as seen in 
the two lowest px bins of Fig. For these bins, no 
excess is observed in the p + p data. In the lowest px 
bin the enhancement in the Au -I- Au data is approxi- 
mately a factor of five above the expectations from the 
cocktail. The data are significantly above the cocktail 
up to ruee = 1 GeV/c^, reaching their maximum around 



TABLE XIL Summary of the fits to the Au -I- Au data 
with the exponential plus the modified power-law function 



{Ae 



-pt/t 



+ B{1 + pj'/b) ") as explained in the text. The 



first and second errors are statistical and systematic, respec- 
tively. 



centrality 


dN/dy(pT > 1 GeV/c) 


r(MeV) 


xVndf 


0-20% 


1.50 ±0.23 ±0.35 


221 ±19 ±19 


4.7/4 


20-40% 


0.65 ±0.08 ±0.15 


217 ± 18 ± 16 


5.0/4 


Min. Bias 


0.49 ±0.05 ±0.11 


233 ± 14 ± 19 


3.2/4 
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rriee ^ 0.4 GeV/c2. 

Figure ^ shows the mass distribution in three px 
bins (0.4-0.6, 0.6-0.8, and 0.8-1.0 GeV/c) in the LMR 
and a possible transition from 1/mee behavior at higher 
Pt (LMR I) to much larger enhancement at lower px 
(LMR II). For the highest pr bin (0.8-LO GeV/c) 
the excess is approximately a constant factor above the 
cocktail. This means that the mass spectrum is still 
close to 1/mee expected for internal conversion. The 
large enhancement seems to appear for the next pT (0.6- 
0.8 GeV/c) bin. For the lowest pr bin the shape appears 
to be different from the 1/mee behavior. 



Figure 36 shows (data-cocktail) //dir(TOee) for the 
three low pT bins. These ratios are proportional to the 
S{mee) factor, and a constant S{mee) leads to a constant 
ratio 7? as a function of mass. While in Fig. ^(a) R is 
still consistent with a constant as a function of mass. 
Figure ^(b) suggests that there is an enhancement for 
0.1 < TOee < 0.4, although the statistical error is too 
large to be conclusive. Figure |3^(c) suggests a large and 
broad enhancement around mee — 0.4 GeV/c^. 

We test whether the R distributions in Fig. ^ are 
consistent with a constant. For each pT bin, we fit a 
constant to the data. The results of the fits are shown 
as the horizont al ba nd in each panel and are summa- 
rized in Tabic pCIII| . For 0.8 < pr < 1.0 GeV/c the 
fit gives good x'^/NDF. Thus the data are consistent 
with the expected 1/mee behavior. The next pT bin, 
0.6 < Pt < 0.8 GeV/c, gives marginally satisfactory 
xVNDF. For 0.4 < pr < 0.6 GeV/c the xV^DF is 
large and the data are statistically inconsistent with a 
constant, suggesting that the electromagnetic spectral 
function is modified at low px- However, due to the large 
uncertainty of the point at mee — 0.4 GeV/c^ the shape 
cannot be well determined. 

The value (R) obtained from the constant fit corre- 
sponds to the direct photon fraction r. The fit value for 
0.8 < PT < 1-0 GeV/c, (R) = 0.177 ±0.032, is consistent 
with the values of r for higher pT shown in Fig. If we 
extrapolate the pr spectrum of direct photons deduced 
from the previous section to lower pT, the expected direct 
photon fraction for pT < 1 GeV/c is ~ 0.17 or less, since 
the spectrum of decay photons is steeper. The R{m) 
values for 0.4 < pT < 0.6 GeV/c are larger than this 
expectation for mee < 0.4 GeV/c^, suggesting that the 
enhancement in the low px region is larger than that ex- 
pected from internal conversion of direct photons. 

In principle, the distribution of R shown in Fig. ^ 
can be extrapolated to mee = to obtain the fraction 
of real direct photons, even if the distribution of R is 
not flat. However, due to large uncertainties in R for 
Pt < 0.8 GeV/c arising from the combination of multi- 
ple dilepton sources, we cannot reliably extrapolate the 
virtual photon yield to mee = to determine the real 
direct photon yield for these two pT bins. 
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FIG. 36: (Color online) Ratio of _R = 

(data — cocktail)//dir('Tiee) for different pr bins ((a): 
0.8 < PT < 1.0 GeV/c , (b): 0.6 < pr < 0.8 GeV/c, (c): 
0.4 < pt < 0.6 GeV/c) in minimum bias Au -I- An collisions. 
The yellow band in each panel shows ztla band of a constant 
fit value to the data points. 



G. pt Spectra for Different Mass bins 



Figure ^ shows the transverse momentum spectra of 
dileptons in different mass windows for Au -I- Au and 



39 



p + p data: 



1 



2ttpt dprdy 



d^N 



2ttpt dpTdydme 



-drrif. 



(42) 



where mi and m2 are the lower and upper hmits 
of the different mass shoes. In the low-mass slices 
(mee < 0.4 GeV/c^) the spectra are truncated at 
low pair-pT due to the single-track acceptance px > 
0.2 GeV/c. The pair-pT cutoff is mass dependent. The 
Au -I- Au spectra have been divided by iVpart/2 in order 
to ease the comparison with the corresponding spectra in 
p + p. The systematic uncertainty due to A'^part (~10%) 
has not been included. In order to avoid the influence 
of e+e^ decays of narrow vector mesons, the mass re- 
gions around the uj meson (0.78 ± 0.030 GeV/c^) and 
the (/) meson (0.1020 ± 0.030 GeV/c^) are excluded. The 
solid curves in Fig. ^represent the expectations from the 
hadronic cocktail, which includes also the charm decay 
contributions. The charm contribution is calculated with 
PYTHIA. Using the random cc correlation makes no differ- 
ence for niee < 0.3 GeV/c^ since the charm contribution 
is negligible at low masses. This difference increases for 
the higher mass bins, leading to spectra lower by ^ 10%, 
20% and 30% for the mass bins 0.3 < rriee < 0.5 GeV/c^, 
0.5 < TOee < 0.75 GeV/c^, 0.81 < m^e < 0.99 GeV/c^ 
respectively, when the random cc correlation is used. 

First we concentrate on the comparison between data 
and the sum of cocktail and charm. In the low-px region 
{pt < 1 GeV/c) all the p + p spectra are consistent with 
the expectations from the cocktail alone for every mass 
window. In the high-pT region however the p + p data 
show a small excess above the cocktail. The Au -I- Au 
data are in agreement with the cocktail in the mass region 
niee < 0.1 GeV/c^. In higher mass bins the Au -I- Au 
data show a large excess b oth a t low- and at high-pT- 

As discussed in section |VE| , we have extracted the 
direct photon yield from the dileptons spectrum in the 
mass range of 0.1 < rriee < 0.3 GeV/c^. The excess in 
this mass range is consistent with internal conversion of 
direct photons. As shown in Figure ^ the direct pho- 
ton component, which appears as a constant R, extends 
to the iriee > 0.3 GeV/c^. Therefore, there should be 
sizable contribution from direct photons in the dilepton 
spectra for niee > 0.3 GeV/c^. The relation between real 
direct photons and virtual photons is presented in Ap- 
pendix Here we use a constant factor S{mee,q) — 1 
to extend the direct photon component to higher mass 



TABLE XIIL 


Summary of a constant 


fit to the ratio data 


shown in Fig. 


1^ The fit range is 0.1 < 


mee < 0.7 GeV/c^. 


PT(GeV/c) 


(R) 


xVndf 


0.8-1.0 


0.177 ±0.032 


7.7/10 


0.6-0.8 


0.198 ±0.033 


16.3/7 


0.4-0.6 


0.293 ±0.040 


21.3/5 



(mee > 0.3 GeV/c^). The dashed curves in Fig. ^show 
the sum of the cocktail, charm and direct photon contri- 
butions to the dilepton spectra for > 1 GeV/c. 

The dashed curves describe the data well for all mass 
bins both in the Au + Au and the p+p data. This in- 
dicates that the excess above the cocktail and charm at 
high {pt > 1 GeV/c) is consistent with the contribu- 
tion from direct photons. It is surprising that the agree- 
ment holds even for mee > 0.5 GeV/c^, where significant 
modifications of the spectral function may be expected 
due to the presence of the vector mesons. However the 
data have large statistical errors for iriee > 0.5 GeV/c^ 
and additional enhancement over the direct photon con- 
tribution is not excluded. The data at high pT are also 
consistent with the cocktail alone for mee > 0.5 GeV/c^. 

In the Au ± Au data, the enhancement over the cock- 
tail is approximately a constant factor for px > I GeV/c. 
It grows towards low-pT- All the Au + Au px spectra 
for every mass bin above 0.3 GeV/c^ seem to indicate 
that the enhancement with respect to the cocktail below 
1 GeV/c is significantly larger than above 1 GeV/c. For 
Pt > 1 GeV/c, the data has a slope similar to the cock- 
tail, as shown by solid curves. For pT < 1 GeV/c, the 
slope of the data is much steeper than the cocktail. 

In order to study this change of the slope in the 
Au + Au data more quantitatively, we subtract the cock- 
tail plus charm from the data and examine the shape of 
the excess. The px spectra are combined in the mass 
range 0.3 < uiee < 0.75 GeV/c^. In this mass range, the 
low-pT cutoff which artificially truncates the px spectra 
at lower mass is avoided. The combined data also have 
increased statistical significance. 

Figure |3^ shows the pair mx—mo spectrum in Au + Au 
for the pair mass range 0.3 < mee < 0.75 GeV/c^. Here 
mx — \/Px + mQ is the transverse mass of the pair and 
mo is the mean value of dN j duiee in the given mass range 
(0.3 < mee < 0.75 GeV/c^ in this case). We plot the data 
as function of mx — mg since invariant differential cross 
sections of hadrons in p+p, p+A, and A+A collisions are 
generally well described by exponential functions in mx- 
Thus the change in the slope can be seen more clearly 
in the mx spectrum. The mx spectrum shows a clear 
change in the slope around 1 GeV/c^. The slope below 
mx — mo < 1 GeV/c^ is much steeper than that above 1 
GeV/c^. In order to characterize the change of the slope 
in the two mx regions, we fit the mx spectrum with the 
sum of two exponentials: 

d'^N , mx , mx 

— = -exp-— ±^2-exp-— (43) 

lirmxdmxdy ii ±2 

where Ai and A2 are the normalization parameters, and 
Ti and T2 arc the inverse slope parameters. 

The result of the fit is shown in Fig. The upper 
solid curve shows the fit function and the dashed and 
dotted lines are the two exponential components. The 
fit gives Ti = 92.0 ± 11.4"*=** ± 8.4'^y"* MeV and T2 = 
258.4±37.3''''^*±9.6">"'' MeV with xVNDF= 4.00/9. The 
two-exponential fit describes the data well. We note that 
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FIG. 37: (Color online) pr spectra of e~ pairs in p + p (left) and An + An (right) collisions for different mass bins, which 
are fully acceptance corrected. Au + An spectra are divided by A^part/2. The solid curves show the expectations from the sum 
of the hadronic decay cocktail and the contribution from charmed mesons. The dashed curves show the sum of the cocktail 
and charmed meson contributions plus the contrib ution from direct photons calculated by converting the photon yield from 



Figure to the e^e pair yield using Eqs. and 



B14. 



the value of T2 , is somewhat higher than, but consistent 
with, the inverse slope of the exponential component of 
the direct photon spectrum, T = 221 ± 19'^'''* ± 19">"'' 
MeV, obtained in the previous section. 

In Fig. |3^ the same Au + Au pT spectrum is fit with 
the sum of direct photons and an exponential function 
in m-r. The exponential function is to characterize the 
low pt component. The direct photon component is ob- 
tained from the direct photon spectrum in Fig. p4, and 
extended to the larger mass region assuming S{m~) = 1. 
We then convert the photon yield to the e+e~ pair yield 
using Equation ^ Thus the direct photon component 
is fixed and the only free parameters of the fits are the 
normalization and the inverse slope T of the exponential 
component. 

The Pt spectrum and the individual components of 
the fit function, i.e. the direct photon component, and 
the exponential component are shown in Fig. B9|. The 
spectrum is well reproduced by the fit and x'^NDF^ 
16.6/11. The systematic uncertainty accounts for the 
uncertainty on the data and the uncertainty (~20%) of 



the cocktail normalization. From the fit, we extract a 
value of T = 86.5 ± 12. 7*^'^* +11.0-28.4"y''' MeV. The 
yield of the low-p^ exponential extracted from the fit 
contributes more than 50% of the yield of the spectrum. 

Both the two exponentials fit shown in Fig. [SSlan d the 
exponential + direct photon fit shown in Fig. |39| show 
that there is a low inverse slope component with T ~ 100 
MeV for rriT — mo < 0.6 GeV/c^. In order to further 
study the mass dependence of the inverse slope, we cal- 
culated the local slopes of the invariant pair rriT spectra 
obtained from the pr spectra shown in Fig. 37, For all 



the mass bins the cocktail and charm are subtracted from 
the data and the average inverse slope (r(mo)) of the ex- 
cess at mass niQ has been numerically calculated as 



(Timo)) 



S(m^ — mo) f (m^rp) Am^j, 



(44) 
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FIG. 38: (Color online) The mr — mo spectrum for the mass 
range 0.3 < rriee < 0.75 GeV/c^ after subtracting contribu- 
tions from cocktail and charm. The spectrum is fully ac- 
ceptance corrected. The systematic error band includes the 
difference in charm yields in this mass range. The spectrum 
is fit to the sum of two exponential functions which are also 
shown separately as the dashed and dotted lines. The solid 
line is the sum. 



where 



1 



2TrmT drriTdy 

1 d^N 
271 pt dpxdy 



(45) 
(46) 



is the invariant spectrum of the electron pairs after cock- 
tail subtraction shown in Fig. ^ 

Figure 40 shows the inverse slopes calculated in two 
ranges, namely < tot < 0.6 GeV/c^ and 0.6 < niT < 
2.5 GeV/c^. For toq < 0.4 GeV/c^ the spectra are trun- 
cated due to the acceptance; therefore, we do not quote 
any slope here. Also for toq < 0.1 GeV/c^ the slope in 
the range 0.6 < tot < 2.5 GeV/c^ is not shown because 
the cocktail subtraction has too large systematic uncer- 
tainty in this region. The solid and dashed lines show the 
inverse slope of the cocktail, calculated in the same way 
as the data, for the same mass ranges. The inverse slope 
parameters obtained from the two-exponential fit as well 
as the exponential + direct photon fit in the same mass 
ranges are also shown. 

Figure ^ suggests qualitatively different behavior for 
low and high tot. In the large tot — toq range (0.6 < 
Tot — Too < 2.5 GeV/c^) the inverse slope (T) is ap- 



FIG. 39: (Color online) The mr — mo spectrum for the mass 
range 0.3 < mee < 0.75 GeV/c^ after subtracting the contri- 
butions from cocktail and charm. The spectrum is fully ac- 
ceptance corrected. The systematic error band includes the 
difference in charm yields in this mass range. The spectrum 
is fit to the sum of direct photons and an exponential func- 
tion. The dashed-dotted line is the direct photon contribu- 
tion. The exponential fit to the low-rriT enhancement is also 
shown (dashed line). The sum is shown with the thick solid 
line. 



proximately 300 MeV and similar to that of the cocktail 
(shown in the dashed-dotted line). In the low tot — toq 
range the inverse slope is approximately 100 MeV, simi- 
lar to that obtained with the two fit methods. This latter 
one (shown at tot — toq ~ 0.45 GeV/c^) is lower than 
that of the cocktail in similar kinematic range (shown in 
the solid line). 

The effective temperature of the lower inverse slope 
component T ~ 100 MeV, obtained from the two fit 
methods as well as the numerical calculation, is much 
lower than the inverse slope of hadrons with similar 
masses (kaons) measured by PHENIX . The slope of 
the kaon spectrum is larger than 200 MeV. The hadron 
slopes rise linearly with mass, consistent with the expec- 
tations from radial expansion of the hadronic source. If 
arising from thermal radiation of the fireball, dominated 
by pion-pion annihilation tt+tt^ ^ p — > e"'"e~, the ex- 
cess yield in the LMR would show similar temperatures 
and a similar linear rise, reminiscent of radial flow of a 
hadronic source |Q. The value of the low-pT inverse 
slope is lower than or similar to the freeze-out tempera- 
ture. Also the inverse slope of dileptons with an average 
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FIG. 40: (Color online) Local inverse slope of the rriT spectra 
of electron pairs, after subtracting the cocktail and the charm 
contribution, for different mass bins. The local slope is cal- 
culated in different mass ranges, < niT — mo < 0.6 GeV/c^ 
and 0.6 < niT — mo < 2.5 GeV/c^. The solid and dashed lines 
show the local slope of the cocktail for the corresponding mass 
ranges. 



mass smaller than 0.5 GeV/c^ is more than a factor 2 
smaller than that of kaons. The simplistic expectation, 
that 0.5 GeV/c^ dilcpton emission is created similarly to 
kaons (from an equilibrated flowing source), is not sup- 
ported by the data. 



VI. THEORY COMPARISON 

The Au + Au e+e^ spectra are now compared to dif- 
ferent models of e+e^ production in the LMR and in the 
IMR. These models, employed at SPS energies, identi- 
fied the pion annihilation process as the main source of 
thermal dileptons in the hadronic phase of the fireball. 
However, this process, mediated by the intermediate p 
meson, failed to describe the observed enhancement in 
the LMR ||8^ at the SPS energy when vacuum proper- 
ties of the p are used. This suggested that in-medium 
modifications of the p spectral function could be respon- 
sible for the enhancement of dilepton yield below the p 
mass. The proposed modifications mostly followed two 
different approaches: 

• The dropping mass scenario followed the scaling 
conjecture of G.E. Brown and M. Rho |Q, which 
postulates that the mass of vector mesons decreases 
in dense matter as the quark condensate (qq) de- 
creases due to partial restoration of chiral symme- 
try. This leads to a decrease of the mass of p meson 
from its vacuum value (0.77 GeV/c^) and causes en- 
hancement of the dilepton yield below the p mass. 



• The broadening mass scenario explains the LMR 
dilepton enhancement by hadronic many-body in- 
teractions The spectral function in a hot 
and strongly interacting hadron resonance gas is 
calculated. The many-body interactions cause 
the broadening of the p resonance, leading to en- 
hancement of dilepton yield below p mass. Other 
hadronic many-body interactions contribute to low 
mass enhancement. 

There seems now to be consensus that the dropping 
mass scenario alone cannot adequately reproduce the 
SPS data. These two scenarios, sometimes in combina- 
tion, are used by different groups to calculate the dilepton 
yield at the top RHfG energy {y/sJjN = 200 GeV). Here 
we compare calculations by the following three groups to 
the PHENIX data. 

• Rapp and van Hees ||l^, |l8, Q calculate the rate 
of dilepton emission from a hadronic gas in thermal 
equilibrium, fn the calculation, the electromag- 
netic spectral function in the vacuum is constrained 
by the data of e~^e~ annihilation into hadrons. 
The spectral function in the medium is modified 
by hadronic many-body interactions. The spectral 
function is characterized by the light vector res- 
onances p(770), w(782) and (/)(1020) at low-mass 
according to the vector dominance model (VDM), 
and a perturbative quark-antiquark continuum at 
higher masses. Dilepton production from a par- 
tonic phase outshines the hadronic gas radiation 
for rriee > 1-5 GeV/c^ at RHIC energy due to high 
initial temperatures. 

• Dusling and Zahed ||l^, ^ use a chiral re- 
duction formalism to calculate the electromagnetic 
current-current correlation function in the medium. 
The experimental data of e+e^ annihilation, r- 
decay, two-photon fusion reactions, and pion ra- 
diative decays are used to constrain the correla- 
tion function. The dilepton emission rates from 
hadronic gas at finite temperature and baryon den- 
sity are then computed from a hydrodynamic evo- 
lution model of Au -I- Au collisions. The partonic 
contribution, which does not become dominant be- 
low the (f) mass, is computed using the Born qq 
annihilation term. 



• Gassing and Bratkovskaya |2^, g7| use a mi- 
croscopic relativistic transport model (HSD) that 
incorporates the relevant off-shell dynamics of the 
vector mesons. This model is well established to de- 
scribe the yields, the rapidity distributions and the 
transverse momentum spectra of hadrons in p + A 
and A + A collisions from SIS to RHIG energies. 
The model reproduces well the dilepton mass spec- 
trum in p + p collisions. In Au -I- Au the model in- 
cludes a modified p spectral function according to a 
collisional broadening scenario as well as a tunable 
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FIG. 41: (Color online) Invariant mass spectra of e^e~ pairs in Min. Bias Au + Au collisions in the IMR. (top left) The 
data are compared to the sum of cocktail+charm. The data are also compared to the sum of cocktail+charm and partonic 
contributions from different models. The calculations are from (center) Rapp and van Hees @ |l|, H and (right) Dusling and 
Zahed jl^, The partonic yields (PY) have been added to the two scenarios for charmed mesons decays, i.e. (i) pythia 

and (ii) random cc correlation. 



dropping p mass scenario. No yield from the par- 
tonic QGP phase is available at the moment from 
the HSD model. 

We received numerical values of these model calcula- 
tions from the authors. The e^e~ rates from these cal- 
culations are filtered into the PHENIX acceptance and 
compared to the data. 

The theory calculation by Rapp and van Hees (l|, |l|, 
|83t is done for a fixed impact parameter 6 = 8 fm. In 
this calculation the number of charged tracks Nch — 230, 
which is consistent with the measured Nch in Au -I- Au 
MB dat a p| . The theory calculation by Dusling and Za- 
hed ^4 I^Sf is done for a fixed impact parameter for 
6 = fm or iVpart = 378. 

Since the calculations were provided for central col- 
lisions, or for collisions with a fixed impact parameter, 
in the comparison to our Min. Bias data we normalize 
the theory calculations by iVpa?f VA^^r't!'^'''", i-e- assum- 
ing that the dielectron yield scales with Nch- However, 
this scaling procedure may have introduced some bias 
in the comparison, as the data show (Fig. ^) that the 



dielectron yield increases faster than A^part (and Nch is 
proportional to A^part)- A more detailed comparison re- 
quires knowledge of centrality dependence of the dielec- 
tron yield both in the data and in the theoretical model. 

For comparison to the data, we add these calculations 
to the hadronic cocktail and the charm decays. The con- 
tribution from the freeze-out p meson is subtracted from 
the cocktail to avoid double-counting. Because the cal- 
culation with the HSD transport model by Gassing and 
Bratkovskaya |2^, ^ can sample any impact pa- 
rameter, the final Min. Bias cross section is obtained 
by performing the integration over impact parameter h 
with a proper geometrical weight. In this case the model 
also calculates the hadronic contributions, which are in a 
good agreement with the cocktail. Only the charm con- 
tribution is taken from the PYTHIA calculation |^ or 
from our calculation with random correlation. 



In Section VIE the Au -I- Au photon spectrum is 
compared to several theoretical predictions. These em- 
ploy hydrodynamical models to calculate thermal pho- 
ton emission from the thermalized partonic and hadronic 
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FIG. 42: (Color online) Invariant mass spectra of e^e'^ pairs in Au + Au collisions in the LMR. The data are compared 
to the sum of cocktail+charm (top left). The data are also compared to the sum of cocktail+charm and hadronic+partonic 
contributions from different models. The calculations are from The calculations are from (top right) Rapp and van Hees jl^ , 
[l^ , ^ , (bottom right) Dusling and Zahed jl^, ^ , and Gassing and Bratkovskaya . 



phases of the reaction, added to NLO pQCD calculations 
that describe prompt photon emission from perturbative 
parton-parton scatterings in the first tenths of fm/c of 
the collision process. 

The experimental conditions reached at mid-rapidity 
in central heavy ion collisions at RHIC of nearly zero net 
baryon density and longitudinally boost-invariance in the 
initial conditions facilitate the applicability of hydrody- 
namics to describe the reaction evolution. In addition, 
the thermalization times usually assumed in the hydro- 
dynamical models {rtherm ^ Ifm/c) are, for the first time 
at RHIC, above the lower limit imposed by the transit 
time of the two colliding nuclei (tq = 2i?/7 « 0.15 fm/c 
for Au + Au at 200 GeV). 

Hydrodynamical approaches describe, under the as- 
sumption of local conservation of energy and momentum, 
the evolution of the system using the equations of motion 
of perfect relativistic hydrodynamics complemented with 
a set of initial conditions (e.g. initial temperature Tinit at 
thermalization time tq), the equation-of-state of the sys- 
tem, and the freeze-out conditions. These models have 
been very successful in describing quantitatively most of 



the differential observables of bulk hadronic production 
(in particular those sensitive to early-times pressure gra- 
dients). 

The same hydrodynamical models, with initial condi- 
tions chosen so as to reproduce the bulk hadron data, 
are now employed to carry out the description of thermal 
photon production over the whole space-time evolution 
of the system. 



A. Comparison in the IMR and constraint on 
possible QGP radiation 

In this Section we compare the model calculations with 
the data in the IMR and investigate whether the experi- 
mental data can constrain the QGP radiation in the IMR. 
We also study whether we can constrain contributions 
from charm and QGP radiation in the LMR. 

All theoretical models predict that there is large contri- 
bution from QGP radiation in the IMR. The QGP radia- 
tion competes with the dileptons from correlated charm, 
which contribute much more than any other cocktail con- 
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tribution in the IMR. There is a large uncertainty in the 
charm contribution in Au + Au, since charm quarks are 
known to suffer energy loss in the medium. Therefore the 
€^ e~ mass shape from semi-leptonic decays of charmed 
quarks niay be modified. The shape calculated from 
PYTHIA provides an upper limit for the expected 
e^e~ yield. The shape with random cc correlation is 
softer for > 0.5 GeV/c^, and this provides the lower 
limit. Below 0.5 GeV/c^ the two shapes are almost iden- 
tical. 

These two scenarios for the open charm contribution 
are added to the predictions for the QGP radiation pro- 
vided by the models described above and are compared 
to the experimental data in Fig. In all the mod- 
els the yield in the QGP phase arises entirely from the 
qq e^e~ annihilation process. The magnitude of 
the yield is closely linked to the thermalization time 
tq. a larger tq translates into a reduction of the ini- 
tial temperature and thus of the total QGP radiation. 
In Fig. ^ the differences between the models are at- 
tributed to different initial conditions used for the hy- 
drodynamic evolution of the QGP phase (ro=0.2 fm/c 
for Dusling and Zahe d |19| , |8^ , |85[ |, =0.6 fm/c or Rapp 
and van Hees |l|, |8||^ 

In the IMR the data have large statistical errors and 
systematic uncertainties. Therefore, they do not allow 
discrimination either between the two proposed scenarios 
for charm production (PYTHiA |0 or random cc corre- 
lation) or among the three theoretical models. 



B. Inclusive Low Mass Excess 

The data in the IMR do allow setting an upper limit on 
the contribution arising from charm or from qq e^e~ 
going to the LMR. Indeed we can saturate the IMR yield 
either with charm or with the partonic yield calculated by 
the theorists and see what their contribution would be in 
the LMR. Thus, neither the charm nor the contribution 
from qq — > e+e" can be solely responsible of the LMR 
enhancement. 

In the LMR the shape and the yields of the e+e" mass 
spectra from charm calculated by pythia are not 
very different from those given by a calculation which as- 
sumes random cc correlation, because the shape is mostly 
determined by the geometrical acceptance. Thus the 
dilepton yield measured in the IMR gives strong con- 
straint on the charm contribution in the LMR. Since 
the calculated charm contribution is consistent with the 
data and is less than the hadronic cocktail below rriee < 
0.5 GeV/c^, we conclude that charm contribution alone 
cannot explain the large enhancement observed in the 
LMR. A similar consideration can be given for the QGP 
radiation. The contribution from qq e^e^ process is 
negligible in LMR I. 

Figure ^ compares the inclusive mass spectrum in 
the LMR with the cocktail-|-charm only and with cock- 
tail-|-charm and the calculations by three groups: Rapp 
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FIG. 43: (Color online) Invariant mass spectra of e"'"e~ pairs 
in Min. Bias Au -I- Au collisions for different pr windows com- 
pared to the exp ectations from the calculations of Rapp and 
van Hees jl^, ^ , separately showing the partonic and the 
hadronic yields and the different scenarios for the p spectral 
function, namely "Hadron Many Body Theory" (HMBT) and 
"Dropping Mass" (DM). The calculations have been added to 
the cocktail of hadronic decays (where the contribution of the 
freeze-out p meson is subtracted) and charmed meson decays 
products. 




and van Hees 1 15 



I , Dusling and Zahed 
and Gassing and Bratkovskaya 
three models the e'^e" yield in the LMR arises mostly 
from the hadronic phase. Rapp and van Hees ^ |8^ 
propose three different scenarios of vector mesons spec- 
tral functions: (i) no medium effects, (ii) dropping p 
mass, and (iii) broadening p mass. Dusling and Za- 
hed 1 19: 1^ ^ use a broadening p mass scenario in 
the hadronic phase and a pion chemical potential of 
p^=50 MeV. Gassing and Bratkovskaya (2^, ||, |8|] 
propose two scenarios: (i) broadening p mass, and (u) 
dropping and broadening mass. 

The common characteristic of the in-medium effects in 
these models is a slight suppression of the yield in the 
p-uj region compared with the unmodified p scenario and 
an enhancement in the region 0.4 < TOge < 0.7 GeV/c^. 
The (f) survives as a pronounced resonance, although its 
width is broadened. These features become less distinct 
once the cocktail contribution and the smooth yield from 
the QGP (which constitutes 15-20%) are added. The dif- 
ferences in the yields of e'^e" pairs in the various models 
are attributed to differences in the medium effects on the 
spectral function, different durations of the lifetime of the 
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FIG. 44: (Color online) Invariant mass spectra of e^e~ pairs 
in Min. Bias An + An collisions for different pT windows 
compared to the expectations from the calculations of Dusling 
and Zahed p^ , p4| , p5[ , separately showing the partonic and 
the hadronic yields. The calculations have been added to the 
cocktail of hadronic decays (where the contribution of the 
freeze-out p meson is subtracted) and charmed meson decays 
products. 



fireball in the hadronic phase, and different evolutions of 
the temperature as a function of time. 

While the calculations proposed by Rapp and 
van Hees m, 



18 



agree with the data for m^e > 
0.5 GeV/c^, the ones of Gassing and Bratkovskaya EG, 



86 



touch the lower end of the systematic uncer- 
tainty in the same mass region. The yields calculated 
by Dusling and Zahed |l^, |5) appear everywhere too 
low to add significant contribution in the LMR, where the 
data are enhanced with respect to the hadronic cocktail. 

All of the models under predict the data for 0.2 < 
niee < 0.5 GeV/c^ by at least a factor of two. It 
should be noted that the contributions in the region 
niee < 0.4 GeV/c^ are very different in the three models 
compared. In Rapp and van Hees jTB], ^ |8^ this con- 
tribution arises from processes like al ttj* — s- Tre^e^ 
or ^ Nj* Ne+e'. Those processes in the HSD 
model are suppressed by a few orders of magnitude and 
are not seen at all compared to the major Dalitz decays. 
In Dusling and Zahed |8^ the main contribution 

below the two-pion threshold comes from 11^, the axial- 
vector contribution in medium. However the absolute 
yield of this process is too low because it is concentrated 
at very low-p^ which is suppressed by our acceptance cut 



FIG. 45: (Color online) Invariant mass spectra of e"'"e~ pairs 
in Min. Bias An + Au collisions for different pr windows col- 
lisions compared to the expectations from the calculations of 
Gassing and Bratkovskaya j20| |86| , separately showing 
the partonic and the hadronic yields calculated with different 
implementations of the p spectral function, namely according 
to coUisional broadening, with or without a dropping mass 
scenario. The calculations which include the dropping mass 
scenario have been added to the cocktail of hadronic decays 
(which is calculated by the HSD model itself) and charmed 
meson decays products. 



(p^"sie > Q 2 GeV/c). 



C. PT Dependence of Low Mass Excess 

Figs. |4^, lU and ^show the e+e^ invariant mass spec- 
tra in different px windows from data compared to the 
sum of all these contributions for the predictions of Rapp 
and van Hees |l|, |l|, ||, Dusling and Zahe d H H pl, 



H, |§r^|J, re- 

the hadronic and 



and Gassing and Bratkovskaya |20 
spectively. The contribution from 
the partonic medium and the charm expectations from 
PYTHIA are shown separately. The charm distribu- 
tion from PYTHIA is somewhat harder than the calcu- 
lations which assume random correlation of the cc pair 
but not very different in the LMR, where the shape of 
the distribution is essentially determined by the detector 
acceptance. 

From the comparison we learn that in general the yield 
from these theoretical predictions is insufficient to ex- 
plain the observed enhanced dilepton production, both at 
low-and high- momenta. At low-px, where the enhance- 



47 



ment reaches approximately a factor of five, the shape 
of the enhancement shown by the data is quite different 
from any of the theoretical models. 

At high-pT- {pt > 1-0 GeV/c) the enhancement is 
about a factor two over the cocktail and its shape is 
quite similar to that of the cocktail. In the previous 
Section we showed that this enhancement can be at- 
tributed to internal conversion of virtual direct photons. 
In the theoretical calculations of direct photon emission 
at RHIC energies, the contribution from the QGP phase, 
e.g. quark-gluon Compton scattering, is the dominant 
source of real thermal photons for px > I GeV/c. The 
process that produces real photons in the QGP should 
also contribute low mass e+e" pairs at high px- How- 
ever, none of the three models includes such processes 
(e.g. q + g^q + j*^q + e'^e~). The QGP radiation in 
these models only include q + q annihilation. This could 
explain the discrepancy between the models and the data 
for PT > 1.0 GeV/c. 



D. Pt Spectra in the Low Mass Region 



the shape of the excess is well described by a pure ex- 
ponential with inverse slope T ~ 220 MeV. If the direct 
photons in Au + Au collisions are of thermal origin, the 
inverse slope T is related to the initial temperature Tjnit 
of the dense matter. In hydrodynamical models, Tinit is 
1.5 to 3 times T due to the space-time evolution Q. 

Figure |4^ compares the direct photon data in cen- 
tral 0-20% Au -I- Au collisions with several theoretical 
calculations of thermal photon emission added to the 
pQCD calculations [|||. Note that the curve by Q 
includes pQGD contributions, while the others do not. 
For Pt < 3 GeV/c, the thermal contribution dominates 
over pQGD. These hydrodynamical models can repro- 
duce the high pT central Au -I- Au data within a factor 
of two. These models assume formation of a hot QGP 
with initial temperature ranging from Tjnit = 300 MeV at 
thermalization time rn = 0.6 fm/c to Tinit = 600 MeV at 
To = 0.15 fm/c lH, 1^ 1^, §|, ll Igl . Figure || summa- 
rizes the Tinit and tq for theoretical calculations shown 
in Fig. 47, There is a clear anti-correlation between Tinit 
and Tq. Lattice QGD predicts a phase transition from 
hadronic phase to quark gluon plasma at ~ 170 MeV. 



The Pt spectra of the excess (i.e. after subtracting 
the hadron cocktail and the charm from the dilepton 
spectra) can be also compared to the theoretical mod- 
els. We already noted that the yields from the par- 
tonic medium in the theoretical models are produced only 
via the qq e+e^ annihilation process. Processes like 
q + g ^ qe^e^ are not included. 

Figure |4^ shows the pT spectrum in the mass window 
0.3 < rriee < 0.75 GeV/c^ after subtracting the contribu- 
tion from the cocktail and the charm. The spectrum is 
compared to the theoretical calculation from Rapp and 



van Hees |18|, |8^ , Dusling and Zahed U 
and Gassing and Bratkovskaya 12^ 



27, R6 37 



B E3' 

respec- 



tively. The figure shows separately the e"*"e~ yields from 
the partonic phase and the hadronic phase (with two 
possible implementations of the p spectral function) and 
their sum is compared to the data. In all the models 



the sum of the cocktail contribution and the 



yield 



from medium-effects is insufficient to explain the exper- 
imental data, and divergences are observed both at low- 
and high-pT- While for Rapp and van Hees ||l^, |l8[ |8^ 
and for Dusling and Zahed [|l9|, HJ, ^ the disagree- 
ment with the data is stron g a t low-p^, for Gassing 
and Bratkovskaya |20, |2^, g7| a better agreement is 



achieved over the full pt range. However, the data seem 
still higher than the theoretical calculations. 



E. 



Theoretical Comparison to Direct Photon 
Measurement 



In section VE we have extracted the direct photon 
yield from the analysis of LMR I. The obtained direct 
photon spectrum in central Au -I- Au, shown in Fig- 
ure 34. shows excess over Taa scaled p + p data, and 



VII. SUMMARY AND CONCLUSIONS 

PHENIX has measured dilepton production in 
Au -I- Au and p + p collisions at ^/snn = 200 GeV. The 
measured e'*"e~ yield is compared with the cocktail of 
known sources of light hadron decays. Gocktail sources 
are mostly measured by PHENIX in the same experimen- 
tal run via hadronic decay channels. Extrapolations to 
low-pT, where experimental data are not always available 
are obtained using an rriT-scahng procedure. 

In the p+p data the e'^e' invariant mass spectrum has 
been measured in the mass range from to 8 GeV/c^ for 
all pair-pT- The intermediate-mass-region (IMR) is dom- 
inated by semi-leptonic decays of heavy fiavor mesons for 
which the extracted production cross sections is consis- 
tent with fixed-order next-to-leading-log (FONLL) pre- 
dictions and with the PHENIX measurement of single 
electrons p^ . 

The low-mass-region (LMR) can be described by 
known contributions from light meson decays, and vir- 
tual direct photons for which the extracted cross sec- 
tion is consistent with NLO pQGD calculations and with 
PHENIX measurements of real photons. 

In Au + Au collisions, the data are consistent with the 
expectations from correlated cc production for Wee > 
0.5 GeV/c^. However, this interpretation is ambigu- 
ous, due to the interplay between possible two differ- 
ent medium effects: energy loss of charm quarks in the 
medium which would deplete the yield in the IMR, and 
QGP radiation, which would increase the yield in the 
IMR. 

In the low mass region the Au + Au Min. Bias inclu- 
sive mass spectrum shows an enhancement by a factor of 
4.7 ± 0.4"*^' ± l.S^y'^' ± 0.9'"°'i<=' compared to the expec- 
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FIG. 46: (Color online) pr spectra of e^e~ pairs for 0.3 < mee < 0.75 GeV/c^ in Min. Bias Au + Au collisions compared 
to the expectations from the calculations of respectively R. Rapp and van Hees jl^, Dusling and Zahed 

Gassing and Bratkovskaya |2^, The spectra are fully acceptance corrected. The curves show separately partonic 

and hadronic yields. For the curves of Rapp and van Hees |l5| , [l^ , ^ the two scenarios: Hadron Many Body Theory (HMBT) 
and Dropping Mass (DM) are shown. The sum is calculated with HMBT. The calculations are compared to the data from 
which the contributions of the cocktail of hadronic decays and charmed meson decays have been subtracted. 



tation from the hadronic cocktail. The enhancement is 
concentrated at low pr {pt < 1 GeV/c). The integrated 
yield increases faster with the centrahty of the collisions 
than the number of participating nucleons. 

At low mass {mee < 0.3 GeV/c^) and high pT (1 
< Pt <5 GeV/c) an enhanced e+e^ pair yield is ob- 
served both in p + p and Au + Au collisions. The mass 
dependence of the excess is consistent with that expected 
for virtual direct photon production. This excess is used 
to infer the yield of real direct photons by extrapolating 
to mee = 0. A perturbative QCD calculation is consis- 
tent with the real direct photon cross section in p + p 
extracted by this method, while in central Au -I- Au 
collisions much larger yields compared with the p + p 
cross section scaled with Taa are observed. In central 
Au -I- Au collisions, the excess over the p + p cross sec- 
tion scaled by Taa is exponential in pr , with inverse slope 
T = 221 ± 19""** ± ig^y"* MeV. 

In Au + Au collisions at very low px there is a further, 
very significant, enhancement that increases strongly 
with centrality. The pr spectrum of dileptons in this 
region has been analyzed with two fit methods and a nu- 
merical calculation. An inverse slope of T ~ 100 MeV 
has been extracted for ttit < 0.6 GeV/c, lower than for 
hadrons with similar mass and similar to the freeze-out 
temperature. 

The Au -I- Au data are compared to different models 
which provide additional e+e~ yield in both the LMR 
and the IMR. In the IMR the data have too large un- 
certainty to discriminate different possible scenarios of 
charm production and QGP radiation. In the LMR no 
quantitative agreement has been found yet with the mod- 
els. 

The yield of direct photons in Au -I- Au collisions is 
compared with several hydrodynamical models of ther- 



mal photon emission at RHIC energies. The models as- 
suming the formation of a hot system with initial tem- 
perature ranging Tinit — 300 600 MeV at times 

To ~ 0.6 0.15 fm/c are in qualitative agreement with 

the data. Lattice QCD predicts a phase transition from 
hadronic phase to quark gluon plasma at T ~ 170 MeV. 

In conclusion, we presented measurements of the e+e~ 
continuum in p+p and Au -I- Au at = 200 GeV in a 

wide range of mass and transverse momenta. In Au -I- Au 
collisions, a large enhancement of the yield of e+e~ pairs 
is observed at low mass and low px in Au + Au. The 
yield of direct photons is deduced from low mass, high 
Pt e'^e~ pairs. Future measurements with an upgraded 
PHENIX detector with higher statistics, together with 
further advance in theory, will allow more detailed study 
of the properties of the hot dense matter formed in heavy 
ion collisions at RHIC. 
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FIG. 47: (Color online) Theoretical calculations of thermal 
photon emission jss], ]9^, [oij are compared with the 

direct photon data in central 0-20% Au + Au collisions shown 
separately and added to pQCD calculations. In contrast to 
the others, the curve by |94] includes pQCD contributions. 
The black solid curve show the pQCD calculation, scaled by 
Taa. The QCD scale fi is set to pr for this calculation. The 
two black dashed curves around the black solid curve show the 
scale uncertainty, with the upper curve and the lower curve 
corresponds to jj. — 1/2 ■ pr and /j. — 2 ■ pr, respectively. 
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FIG. 48: (Color online) "Unit vs. tq for various theoretical 
calculations shown in Fig. 
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APPENDIX A: BACKGROUND 
NORMALIZATION 

1. Pairing of Electrons and Positrons 

In the following we assume that, as dictated by the 
charge conservation law, e~ and e"^ are always produced 
in pairs and that most of these pairs are produced statis- 
tically independent of each other. Let us say N pairs are 
produced in a particular event and N is given by a prob- 
ability distribution P{N). Of the N pairs only a fraction 
Ep is fully reconstructed, and then the number of recon- 
structed pairs Hp is given by a binomial distribution B 
sampling out of N "events" with a probability Sp. 

• Probability to get Up pairs from N true pairs: 
uj{np) = B{np,N,ep) 



with an average: 



2p) = EpN 



• and variance: cr^ = epN{l — Sp) 

Of the remaining pairs one track is reconstructed with 
a probability e_|_ or e_. For a given N and Up the num- 
ber of additional single positive tracks n+ and negative 
tracks n_ follow a multinomial distribution M with three 
possible outcomes for each of the N — Up unreconstructed 
pairs: no track, one -I- track or one — track. 

The probability to get n+ and n_ single tracks from N 
true pairs with Up reconstructed pairs, i.e., from (N — Hp) 
not fully reconstructed pairs is: 

a;(n+, n_) ~ M{n^,n^; N — Up, e+, e_) 

N-Up 

uj{n+) = ^ M{n+,n^;N ~ np,e+,e-) 



M{n+,n^;N ~ np,e+,S-) (Al) 



• with average: {n±) ~ e±{N — Up) 

• variance: cr^ = e±{N — np){\ — e±) 
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and covariance: cov(n+,n_) = — (N — np)£+£_ 



In this case the number of unUke-sign pairs for a given 
and Up is: 

N—rip N — Up 

(n+_) = nl + rip ^ n+a;(n+) + np ^ n_a;(n_) 

n+— 1 n_=l 

N — rip N—Up 

+ X] X] n+n-i^{n+,n-) 

n+ = l n_ = l 

= rip + npS+{N - rip) + npe-{N - rip) + (n+n_) 
= rip + npe4.(iV - rip) + npe_(iV - rip) + 

£+£_ (iV - rtp)^ - £+£_ {N - rip) 
= n,p + e+Nup - e+Up + e-Nup - S-Tip 

+£+£_ A^^ — 2e+e-Nnp + e+e-Up 
—e+E-N + £+£_np 
= (rip + £+ (A'' - rip)) (rip + £_ (AT - Up)) - 

s+S-{N-np) (A2) 

Similarly we can calculate the number of like-sign 
pairs: 

N-rip 

2(n++) = (rip + n+)(np + n+ - l)a;(n+) 

n-i- — 1 

= rip - rip + (n^) - (n+) + 2np(n+) 

= - rip + £^(Af - np)2 + £+(1 - £+)(A- - rip) 



-£+(A' - rip) + 2£+np(Af - rip) 
= nl-np + e\{N - Upf - e\{N - rip) 
+2s+np{N - Up) 



(A3) 



and 



2(n__) = nl - Up + e'i{N - Upf 

_2 



-£i(Ar-np) + 2£_np(A''-np) (A4) 

To obtain the expected number of like- and unlike- 
sign pairs for a fixed number of real pairs A'' we need to 
average over all possible reconstructed pairs Up-. 

(iV+_) = J2i^+_)B{np,N,ep) 

Up 

+{e+N + e-N - 2s+S-N + s+e-){np) 
+e+e-N^ - e+e-N 

= (l-£+-£_-h£+£-)(£pAr^-t-£p(l-£p)iV) 

+{e+N + e-N - 2e+S-N + s+S-)epN 
+e+e-N'^ - e+e-N 

/ 2 2 2 I 2 I I 

— (,£p ~ ^p^+ ~ ^p^— ' ^p^+^— + £p£+ + £'p£ — 

-2£p£+£_ + £+£_)(A^2 - A^) + £pA^ 

= (£p + £+(l-£p))(£p+£_(l-£p)) 

(Af 2 -N)+ epN (A5) 



Now we calculate the like-sign background: 

2(Ar++) = Y.2{n++)B{np,N,ep) 

Up 

= elN^ + £p(l - ep)N - SpN + ele\N^ 
+£+ep(l - ep)N - 2e\epN'^ + e\N'^ 
-e\N + £^£pAf + 2e+£pN'^ 
-2e+elN^ - 2e+ep{l - ep)N 

= el{N^ -N)+ eleliN^ - N) + elepN 
-2e\epN'^ + e\{N'^ -N)+ e\epN 
+2e+epN'^ - 2e+e%N'^ - 2e+epN + 2e+elN 

= [el +el+ elel){N^ - N) - 2elep{N^ - N) 

+2£+£p(Af2 _ _ 2£+£2(Ar2 _ N) 

(Ar++) = l(ep + £+(l-£p))2(Ar2_7V) (A6) 
and 

{N-) = i(£p + £_(l-£p))2(iv2-7V) (A7) 

Finally wc need to average over all N to get the fore- 
ground unlike-sign pairs: 



{FG+-) = J2iN+-)PiN) 



N 

= (£p + £+(l-£p))(£p+£_(l-£p)) 

.i{N^)-{N))+ep{N) 
= {BG+-) + {S) (A8) 

The unlike-sign foreground FG_| consists of the sum 

of the unlike-sign background BG+- and the signal S = 
ep{N). Similarly the like-sign foreground is calculated 
as: 



N 



{FG-) 



-{ep + e+{l-ep)f{{N^)-{N)) 

{BG++) 

Y,{N-)P{N) 



(A9) 



N 



= {sp+e-{l-ep)n{N^)-{N)) 



^ {BG-) 



(AID) 



The like-sign foreground contains no signal. 

So due to the fact that electrons and positrons are al- 
ways created in pairs, the unlike-sign baekgrouncl is the 
geometric mean of the like-sign backgrounds, indepen- 
dent of the primary multiplicity distribution. 

{BG+-) = 2^{BG++){BG—} (All) 
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Let us compare the background to the product of the 
average track multiphcities. For a fixed Up-. 



= (np + n+)uj{n+) 

n^ — l 

= np + (n+) 

= Up + e+{N - Up) 

averaged over all possible Up-. 

N 
np=0 

= EpN + e+N - e+EpN 
= {ep+e+{l~ep))N 

or averaged over all possible N: 

{FG+) = J2{N+)P{N) 

N 

= {ep + e+il-Sp)){N) 

and thus: 



(A12) 



(A13) 



(A14) 



{FG+){FG-) = {ep + e+(l - £p))(£p + £_(1 - ep)){Nf 

(A15) 

or 



{FG+){FG^) 



= 1 



(A16) 



So in general {BG+-) ^ (FG+)(FG_), except for the 
special case that P{N) is a Poisson distribution. Note 
this is the opposite conclusion one derives in the case 
that the sources of + and — tracks are independent, i.e., 
+ and — tracks are produced as singles and not as pairs 
as in the case of muons. In that case {FG+) (FG_) is the 
correct background normalization. 



APPENDIX B: RELATION BETWEEN REAL 
PHOTONS, VIRTUAL PHOTONS, AND 
ELECTRON PAIRS 

1. Introduction 



Figure 49 illustrates that, in general, any source of high 



energy real photons can also emit virtual photons which 
materialize into electron pairs. On the left side a real 
photon is emitted by a source labeled as M{Q^ = 0). On 
the right side is an analogous diagram, where a virtual 
photon with mass m^. is emitted. The virtual photon can 
then convert to an e+e— pair if m-y* > 2me. This e~^e~ 
pair production process is a QED correction to the real 
photon production process and is often called internal 
conversion. 

In the energy region, where electroweak effects are neg- 
ligible, an electron pair can only be produced through a 



M(Q^=0) 




M(Q2=m^) 



FIG. 49: Diagram for real photon production (left) and its 
associated process producing an e^e~ pair (right). 



virtual photon. (Here we don't include e+e^ pairs from 
correlated weak decays such as cc — > e~^e~.) Thus any 
electron pair production process can be described as pro- 
duction of a virtual photon and its subsequent decay into 
an e"'"e~pair. 

The yields of virtual photons dN^. and electron pairs 
dNpp are related: 



dAP 
L{M) 



a L{M) 
3^ M2 



dN. 



7* ' 



Ami 
M2" 



2mi 
M2 



(Bl) 
(B2) 



where M is the mass of the virtual photon or the elec- 
tron pair [M — m.^. — rriee) and a is the fine structure 
constant (a ~ 1/137). The factor, ^3^, is a univer- 
sal factor describing the decay of the virtual photon into 
an e+e~ pair. This relation is exact to first order in the 
electromagn etic coupling a. 
Equation Bl can be written as 



d^N,, 



a L{M) 
'MP ~ 3^ M2 



S{M, q)dN^. 



(B3) 



Here we have introduced S{M,q) = dN^*{M) / dN^ to 
factor out the difference between real photon emission 
and virtual photon emission. The factor S{M, q) is pro- 
cess dependent and accounts for effects such as form fac- 
tors, phase space, and spectral functions. S{M,q) ap- 
proaches 1 for smaU M, S{M,q) 1 for M ~* 0. Ad- 
ditionally, since L{M) ~ 1 — Qm^/M'^ for <^ M, 
L{M) = 1 is a very good approximation. Thus the re- 
lationship between the electron pair yield and the direct 
photon yield simplifies to 



dM^ 

d^N,. 



a 1 

TrdN^, 



dM 



2a 1 
3^M' 



dN^ 



(B4) 
(B5) 



The relation between real photon production and elec- 
tron pair production shown by Equation B5 is valid if 
M <ti Fry, i.e. if the virtual photon is quasi-real. In 
this region, the yield of electron pairs in the mass range 
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mi < M < 1712 is related to the photon yield as 



cross section for virtual photon production as follows p& 



— — dTV^ (B6) 



l^log^d^,. (B7) 



In the following, we discus two examples of internal 
conversions: Dalitz decays and high pT Drell-Yan pro- 
duction. We then discuss the relationship between direct 
photons and electron pairs from thermal sources. At the 
end, we illustrate the relationship between electron pairs 
and virtual photons using a theoretical model calculation. 



2. Dalitz decays 

Dalitz decays of pseudo-scalar and vector mesons are 
prime examples of internal conversion. In these pro- 
cesses, a virtual photon, instead of a real photon, is emit- 
ted in the decay of a hadron and subsequently decays 
into an e^e~ pair. The relation between A — > Be^e~ 
and A is given by Q 



dip 



a L{M) 
Skb{M) = |FAB(Af2)|2x 



Sab{M)V{A^ B-f), 



(B8) 



(1 



-f 



(m 



3/2 



where ttia and tub are the mass of hadrons A and B, 
and Fab is the electromagnetic transition form factor. 
Sab{M) here is an example of S{AI,q) in Equation B3 



For decays of pseudo-scalar mesons (P = n^,ri,r]') the 
relationship between the photonic decay {P — > 77) and 
the corresponding Dalitz decay (P — > e'^e~"f) is given by 
the well known KroU- Wada formula |f73, |74| . 



dr(P-^e+e-7) _ 2a L{M) 
~ 3^ Af2 



SKwiM) 



5Kw(Af)r(P^77), 



ap 

* Tod 



(BIO) 



where mp is the meson mass and Fp{A'P) is the elec- 
tromagnetic form factor. Note that the factor 2 in |^ 
accounts for the fact that each of the two decay pho- 
tons can convert to an electron pair. The form factor is 
usually parameterized as Fp{Q^) = 1/(1 — Q^/Ap). Ex- 
perimental measurements of the transition form factor by 
Lepton-G j?^ and Cello j95| show Ap ~ Mp, consistent 
with the vector meson dominance model (VDM). 



d^c 



pp- 



*i+i- 



X a L{M) 



dM^dp^dy Btt M2 



dp^dy 



(Bll) 



Where M, px, and y are the mass, the transverse mo- 
mentum, and the rapidity of the virtual photon. For 
Pt 3> A/, the virtual photon cross section becomes 
equal to the real photon cross section {da-f* — > da^ as 
M/pt 0). Thus the electron pair cross section can be 
described as 



da„ 



-X 



a L{M) da< 



pp-^jX 



dM^dp^dy 3tt M2 dpj,dy 



(B12) 



Direct photon production via gluon-Compton scatter- 
ing {q+g — > 9-1-7) has an associated electron pair produc- 
tion process {q + g ^ q + ^* q + e~^e^). For the lowest 
order pQCD calculation, the following relation between 
the two processes is obtained by an explicit calculation. 



(i (Tee a L{M) da^ 
dM^dt " 3^ A/2 Hr " 



xSggiu,t,s), (B13) 



where 



Sqg{u,t, s) = (1 + 



2u 



f2 + s2 



A/2) 
2x 



{x + Vl + 2:2)(3a;2 + 1 + 2xVl + a;2) 



(.B14) 



(B9) Here x = pr/M and s, t, u are the Mandelstam variables 
defined as s = {p + fc)2, t = {p — fc')2, u ~ {p ~ p'Y 



p, fc, k' ,p' are 4-momenta of the incoming quark, the in- 
coming gluon, the outgoing (virtual or real) photon, and 
the outgoing quark, respectively. The factor Sqg{u,t,s) 
accounts for the difference between the virtual photon 
cross section and the real photon cross section, and 
Sqg{u, t, s) becomes unit y for small A/. It is an example of 



5(A/, q) in Equation B3. For 90° scattering in the c.m.s.. 



Sqg ~ (l-p|,/5A/2). Thus Sgg ^ 1 US p^ / 0. This 
means that the approximation Sqg ~ 1 is valid as long as 
p'^/M^ <C 1 even if A/ is relatively large. 

In a kinematic region, where pt/M is not very large, 
contributions from parton fragmentation into direct pho- 
ton (real and virtual) become significant. Unfortunately 
hadronic effects in parton fragmentation into photons 
are poorly understood. Kang, Qui and Vogelsang dis- 
cussed Pm the theoretical uncertainties in this process. 
However, the uncertainty due to this effect is relatively 
small, except for very low ij>T ^ 1 GeV/c) [M. 



Thermal radiation 



3. High PT Drell-Yan process 

Irvp + p collisions, the cross section for Drell-Yan elec- 
tron pair production can be expressed in terms of the 



In heavy ion collisions, thermal radiation from the hot 
and dense matter formed in the collision may contribute 
to both real direct photon production and electron pair 
production. The emission rate of electron pairs per space- 
time volume from a thermal source can be described 
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in terms of the electromagnetic (EM) spectral function 
as i6|,|97| 



dRe 



L{M) 



d^q 37r3 M2 

L{M) 
^ ^37r3 M2 



Imnj^„^^(A/,g;T)/^(go,T)(B15) 

im(2nL + nD/''(go,r).(Bi6) 



Here 11^^ and 11^^ are the transverse and the longitu- 
dinal components of the in-medium photon self-energy 
tensor n^„^^, and f^{qo,T) = l/(e«°/'^-l) is the Boltz- 
mann factor. 

Using the same notation, the emission rate of virtual 
photons is described as 



dRy" 



-^Imn^f^,^(M,g;T)/^((Zo,T). (B17) 



The real photon emission rate is obtained in the limit 
of M — > 0. The longitudinal polarization contribu- 
tion Ilg^ vanishes for real photons: n^„(M, q;T) 
n^„(0,g;T),n^„(M,g;r) ^ for M 0. The real 
photon emission rate is given by 



qo 



dRy 



-^Imni,„(Af = 0,q;T)f^{qo,T). (B18) 



The virtual photon and electron pair emission rates are 
thus related by 



dR, 



qo 



dA'Pd^q 



1 dR, 
2 



dR, 



d^q 
a L{M) 
3^ AP 
a L{M) 



dR. 



qo 



7* 



d^q 
dR^ 



90- 



dAPd^q ~ 3n AP d^q 



(B19) 
(B20) 



Note that Equation B19, describing the relationship be- 
tween the virtual photon and electron pair emission, is 
exact to order a in QED and is exact to all orders of the 
strong coupling constant. Note also that this relation- 
ship between electron pairs and photons is the same as 
that shown in Equation Bll , This reflects the fact that 



an e~^e pair can only be produced through virtual pho- 
ton and that the conversion rate of a virtual photon into 
an e^e^ pair is described by a universal factor, ^ ^m^^ • 
Equation B2C is equivalent to Equation B12 and is an 



approximation for small AI, where dR* ~ dR^. 

The relations above are for the emission rates per 
space-time volume. The yields dN^e and dN^ are ob- 
tained from space-time integral over the rates. 



dAP 
d^Nee 

dM 



2a L{A1) 



37r M 
2a L{AI) 
3^ M 



-dN^ 



S{M, q)dNy 



(B21) 
(B22) 
(B23) 



The se a re the same equations as Equ ation Bl - Equa- 



tion P33| . Here, following Equation |B3| , the difference be- 
tween real and virtual photons is factored out in <S'(M, q). 
The S{AI, q) factor for radiation from thermal sources 
can be written as 



S{M,q) = 



(Im(2HL(A/, q) + nL(M, g))/^(<Zo)i 



(I: 



m2Hj„.(0,<7)/S(<Zo;M = 0)) 



{B24) 



Here () indicates the space-time average. 

Deviation of S{AI, q) from unity can arise from 
non-zero Hg„j {AI, q) in the medium and a change of 
H^,„(M, 5) from H^,„(0, g). The behavior of H^,„ and 
Hg^ is model-dependent. However, on very general 
grounds we can conclude that the in-medium spectral 
function Hem is a smooth function of AI for the low-mass 
region {AI < a few 100 MeV/P). 

Hadronic interactions yield and almost flat behavior 
in S{AI) as we see in a model calculation by Rapp later. 
q+g scattering gives an almost flat contribution in S{M); 
see Sqg of gluon Compton scattering. As we see later, 
the contribution from qq a nnihil ation is not constant 
and is cx (see Equation B27). This means that it 
is strongly suppressed in the low-mass region relative to 
hadronic and qg scattering contributions. Thus S{AI,q) 
is a smooth and almost constant function of M at low 
masses. Furthermore, the short time (< 1 fm/c) between 
rescattering among hadrons and partons in the medium 
should smear any feature in S{M) smaller than a mass 
scale of a few 100 MeV/P. 



5. Conversion of e+e pairs to virtual photons 



Equation (B21) can be used to determine the virtual 
photon yield from the electron pair yield. The equation 
can be rewritten as 



dN^ 



90- 



S{A4, q)qo^ 

37r dNee 



(B25) 
(B26) 



Using this relationship, the virtual photon spectrum as a 
function of mass AI and pt can be determined from the 
double differential electron pair spectrum. In particular, 
the mass dependence for a given pT bin can be mea- 
sured. This is a direct measurement of the shape of the 
S{AI, q) function and the shape of (ImH^™(M, T)f^{T)) 
as a function of AI . The real photon yield is then ob- 
tained by extrapolation of dN^* {AI) to the M ^ limit. 



6. Electron pairs and virtual photons in a 
theoretical model calculation 

Figure ^ illustrates the relationship between the elec- 
tron pair mass spectrum and the virtual photon cross 
section. Figure ^ (a) shows the double differential 
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FIG. 50: (Color online) (a) Electron pair emission rate and (b) virtual photon emission rate calculations at a fixed pair 
Pt = 1.525 GeV/c jisj. The solid curve shows the hadronic many-body theory in the medium. The dashed curve shows the 
calculation when the EM spectral function in the vacuum is used. The dotted curve shows the qq annihilation contribution. 



electron pair spectrum, {l/pT)dNee/dMdpTdy at pt = 
1.525 GeV/c from a model calculation of electron pair 
production by Rapp |83|. The dashed and soUd curve 
show electron pairs from the hadronic gas. The calcula- 
tion shown in the dashed curve uses the spectral func- 
tion Hem that is unchanged from its vacuum value so 
the line shapes of vector mesons {p,Lu,(l)) are unmodi- 
fied. The calculation shown in the solid curve, uses a 
spectral function calculated by a model Lagrangian of 
hadronic many-body interactions, and the line shapes of 
vector mesons are broadened due to the interactions. It 
also includes the contributions like ai(1260)^ tt + e+e^, 
p ^ TT + e+e", and -|- tt ^ A^* ^ iVe+e" that fill 
out the low-mass regions below two-pion threshold. In 
the low-mass region, the mass spectrum steeply increases 
with decreasing M. This steep behavior is due to the 
1/M factor in 7* — > e+e^. The dotted curve shows the 
contribution from the leading order (LO) qq annihilation 
which is augmented with a g = Hard Thermal Loop 
(HTL) correction factor ||9|. 

Figure |5^(b) shows the same calculations presented as 
the differential yield of virtual photons as a function of 
mass. The electron pair yield shown in Fig ^(a) is con- 
verted to the virtual photon yield using Equation B26. 
In this plot, the solid curve becomes almost constant for 
M < 0.3 GeV. The steep 1/M behavior of the electron 
pair spectrum is removed, and much more smooth be- 
havior of the virtual photon spectrum is revealed. The 
plot shows that the virtual photon yield is almost con- 
stant as function of AI. The value of the solid curve at 
M = corresponds to the real photon yield. Thus in 
this model calculation, the yield of virtual photons for 
0.1 < M < 0.3 GeV is almost the same as that of real 
photons. The solid curve illustrates that in a consistent 
theory calculation the yield of virtual photons is a smooth 
function of M and it becomes the real photon yield in the 



limit of M = 0, dN*{M,pT) dN^ipr) as M 0. 

The flat behavior of the solid curve in this low-mass 
region comes from the fact that hadronic scattering pro- 
cesses such as ai(1260)-^ n + e+e^, p ^ tt + e+e^, 
and -I- TT ^ N* — > Ne^e^ dominate this low-mass re- 
gion. These processes are internal conversion of the cor- 
responding real photon production processes, i.e., oi — > 
7r-t-7, p— >7r-|-7, 7V-|-7r^A'^-|-7. Thus production of 
virtual photons from these processes should be very close 
to that of real photons at low mass. In the language of 
the spectral function Hem, these processes contribute to 
Hem at M = as well as M > 0. Their contribution to 
Hem is a smooth, almost a constant function of virtual 
photon mass M. 

In Fig. pO|(b), the contribution from qq annihilation is 
shown by the dotted curve. In perturbation theory, the 
qq contribution is given by 

ImHemiM) = ^(1 + ^ + ...)iV,E(e,)2. (B27) 

The quark annihilation contribution behaves as cx 
in the low-mass region. Thus it is strongly suppressed 
and has little contribution in the low-mass region. In the 
high-mass region, the Af^ behavior of the quark annihi- 
lation is suppressed by the Boltzmann factor. 

There is a large uncertainty in the approximation used 
in the dotted curve, and it is shown here just to illus- 
trate that the qq annihilation contribution is sub-leading 
contribution in the low-mass region {M < a few 100 
MeV/c^). Many effects can alter the shape and mag- 
nitude of the qq contribution. The calculation shown 
in Fig. ^ uses zero quark mass. The effective quark 
mass in the medium is uncertain, but in most theoreti- 
cal calculations it is the order of the temperature of the 
medium. Quarks in the medium should also have a large 
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width (more than a few 100 MeV/c^) due to the short 
time between interactions. A non-zero effective quark 
mass would further suppress the quark annihilation con- 
tribution in the very low-mass region while a finite quark 
width would smear the behavior and can cause non- 
zero qq contribution at M = 0. 

It should be noted that the dotted curve does not 
include processes like q + g — > 9 + 7* that are associ- 
ated with real direct photon production. This is be- 
cause HTL calculation of thermal radiation from QGP 
is only available in the real photon case. Contributions 
from processes associated with real photon production 
can be much larger than those from the LO qq annihila- 
tion shown in Fig. |5^ in the low-mass region. Turbide, 



Gale, and Rapp |Q calculated real photon production 
in a hadronic gas using the same model and compared 
it with real photon production in the QGP phase using 
the complete leading order HTL analysis. They found 
that real photons from the QGP outshine those from a 
hadronic gas for > 1-5 GeV/c in Au + Au collisions 
at RHIG. Since the virtual photon yield should be con- 
tinuous from M — (i.e. real photons) to M > 0, this 
implies that the contribution from the QGP, including 
processes like q + g q + j*, is as large as or even 
greater than that shown in the solid curve in Fig. pO|(b) 
in the low-mass and high-pT region (i.e. M < a few 100 
MeV/c^ and pr > 1.5 GeV/c^) at RHIG energies. 
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